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1. INTRODUCTIDH 


A radiative heating and cooling algorithiu has been developed for 
use In the dynamical model of the middle atmosphere written by J. R. Holton 
{Holton and Wehrbein, 1979; Holton and Wehrbein, 1980) • A Curtis matrix 
Is used to compute cooling by the 15p and lOp bands of carbon dioxide. 

Escape of radiation to space and exchange with the lower boundary ate 
used for the 9.6p band of o^one, Voigt line shape, vibraticnal relaxation, 
line overlap, and the temperature dependence of line strength distributions 
and transmission functions are incorporated into the Curtis matrices. 
Properties of this algorithm have been outlined in Wehrbein and Leovy (1981) 
This report is a more detailed description of its development. 

Section 2 discusses the distributions of the atmospheric constituents 
included in this algoritnm. Section 3 describes the method used to compute 
the Curtis matrices. Cooling or heating by the 9.6y band of ozone ts 
discussed in section 4, A description of the FORTRAN programs and sub- 
routines that have been developed in this study is given in section 5. 
Listings of these routines are found in the appendix. 
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2 . ATMOSPHERIC CONSTITUENTS 

With the exception of molecular oxygen p which absorbs solar radiation 
In the upper mesosphere, the atmospheric constituents that determine the 
thermal structure of the middle atmosphere are exceedingly minor constituents - 
A discussion of the contributions of each gas la found in London (1980), In 
the present work we shall consider only carbon dioxide (15y and lOy bands), 
and ozon**, (9,6y band and solar absorption), and molecular oxygen (solar 
absorption), A more sophlstiGatC:d algorithm might include the overlap of 
carbon dioxide and ozone bands with each other, and with the bands of 
tropospheric water vapor when computing radiative exchange with the ground 
(e,g,p Ramanathan, 1976), 

2.1 Carbon dioxide 

Carbon dioxide is, for our purposi?,s, mixed throughout the atmosphere 
with volume mixing ratio of 0.33 x 10""^. Therefore the amount of CO2 between 
any two pressure levels la always the same, and the transmission between two 
pressure levels depends only on the temperature of the intervening material. 

2 . 2 Ozone 

Both the local density and the column abundances of ozone and moxecular 
oxygen are required to compute the radiative heating due to absorption of 
solar radiation. The specific ozone density ~ ’^STP * 

2.6869 X 10^^ and is the number density of ozone. The mass mixing 

ratio (in grams/gram) is given by mass 

densities of ozone and air, respectively. The specific density is related to 
the mass mixing ratio by 


3 


^^3 



V “3”STP 


( 2 . 1 ) 


where p and T are local pressure and temperature, R is the gas constant for 

ft 

air, and a., la the mass of a single ozone molecule. 

.o' 

The specific column abundance (in cm at STP) Is given as 


^3 “ I 

z 


1 

®3"STP^ 


P 

0 


W^(p')dp’ 


( 2 . 2 ) 


Above a certain altitude z, (pressure level p, ) we assume that ozone la 

A 1 

distributed with a constant scale height so that 


njCzj)- - 


STP J 


n(zj)e 


^ (z — z j^) /Hj 




Belov p, 


the meftn value of W« 


J-1 




in each layer is used to obtain 

IWjCPj) +W3<Pj^^)J(Pj^^ - Pj) 


For arbitrary p is interpolated linearly between levels: 

P *■ Pj 

H 3 (p)- n 3 (Pj)r (1 - Y)n 3 (Pj^j_) where y => (2.3) 

J"^l JT 

for Pj i P < Pj,,j . 
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2.3 Molecular oacygen 

Assume that the mixing ratio of oxygen In air is 0.2Q93 up to 155 km, 
(3.645 X, 10”* atmospherr«s) , above which oxygen follows the profile 

n 2 C*) ■ 1 * 10" e ’ . 

T where is Boltzmann's constant. 

The column abundance below 155 km is therefore given ' y 

z 

Do * ! n_(z')dz' - 1.711 X 10^® + (p - 3.645 x 10~* p ) (2.4) 

Liz in 8 o 

0 ^ 

where m Is the mean mass of a molecule of air, p is standard surface 
■ a ■ o 

pressure, and is given in molecules cm“^. 


Below 155 km a^(z) - 0,2095 p/k^ 
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3. DERIVATION OF GURTIl? MATRIX ELEMENTS 

Th€ net upward flua of thermal radiation at wavenumber V Is given 
by 

F^(t) - 2ff f du |j^(t)e“l^ d|t -Ti 

O 0 


- 2it 


,1 
dy 
o o 


fj’(t)e"l‘'’'^i^^ djt-Tl , 

(Goody, 1964, p. 49 ) 


(3.1) 


where |j; Is the absolute value of the cosine of the zenith angle, and the 
optical depth is defined as 


T - t(p) 



k^(p”)o(p”)dp” , 


(3.2) 


where is the absorption coefficient and c the concentration of the 
absorbing gas. The source function is assumed to be isotropic in the 
upper (-♦-) and lower (-) hemispheres. Changing the variable of integration 
to pressure yields 


i .P 




- 2ti 


I dp j J*(p) «xp I - I 


k £S£lll k £<ip. 

V gu / V g P 


0 p 

1 ® 


+ 2tt 


P’ 


dp 
0 P 


. <<>’.) {- 1 " 


k ^dp' 
V 8P i P 8 


s 

1 0 


- 2tr 


dp 
0 P 


j;(p') exp j k^ I dp’, (3.3) 
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where 1b the surface pressure. The second terra, arising from the surfacej 
can be Integrated over p' directly. 

Consider the flux of radiation integrated over a frequency interval Av 
sufficiently snsall that JJ* can be approximated by mean values over frequency 
J*. Then changing the order of Integration gives 


dvF^(p) 

Av 


2ir 


dp'j"^ Cp') 


dn J 
'0 Av 


dV exp 


P’ 

P 


V 


cdlL!4 
gU / 



+• 2rrJ^(p ) 
*^s 


dy y 


d\> exp 


AV 



£d£!4 

gy / 


- 2tt 


dp’j'(p’) 


dy 


exp 


Av 


cdp 
^ gy 




V g 


(3.4) 


Note that each of the integrals in the exponents are non-negative. 
Define the flux equivalent width 


where 


e(p,p’) 


dy y 


dvu 


T^(p,p' ,y)] 


Av 


T^(p,p',y) 



k c/gy dp 

V 



(3.5) 


Function e(p,p') Is simply one minus the transmission function between p and 
p' weighted by y Integrated over one hemisphere. Clearly e(p,p) =» 0, and 
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3P 
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j I «p|- IJ 




t • 
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^ r p' > p 

. (3.6) 
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)r i 

Wi p' < p 


0 Av P 

Changing vatlables froia p to a» th« expression for the integrated flux 
becomes 


dvF^(p) 


Av 


e(P»P_) 

2t j 
0 


/(p*)de(p,p') + 2 tt/(p^) 
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eCPkPg) 


0 


+ 2tt 


j“(p’)de(p,p') 


(3.7) 


e(p,0) 

Define the upward flux 


F (p) ” 2 tt 


e(p,p„) 

s 


j"^(p*)de(p,p') + 2 tiJ'^(p ) 


T~ ^<P'Ps^ 


and the downward flux 


e(p,0) 


f”(p) - 2tt 


j”(p’)de(p,p') 


(3.8) 


Assume the source function is Isotropic over ail directions J = J > and 
define a nondlmenslonal "temperature" 0 =* J/B, where B is the Planck function 
in Av evaluated at some reference temperature. We have 

e(P»P„) 

I dvF^(p) - F^(p) - f”(p) - 2ttB 


6(p’) de(p,p*) 


Av 


+ 2nB 

Av . 

[2 ■ 

- 2ttB 



0 


0(p')de(p,p*) . (3.9) 
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3,X Vertical 

Divide the etaosflt'isre into N layers (see 
is zero, and «t the ground* Pressure 
levels within layers are denoted with catets. 


Fig, 1), The pressure at 
and temperature evaluated at 
To simplify notation, let 


e(p^,Pj) ■ .y 


e(Pi»P|) 




e(Pi,P^) 


®i,N+l • 


For the simplest case 0 is treated as a constant throughout each layer. Then 


N 

^i * S (*ipj+l ■ (t " ®i,N+l) Vl 

j-1 

N+1 

” 2TrB ^ (*i,j+l ■ ®ij)®j where ~ ’ (3,10) 

j-i 

and 

x-1 

*’i * - ^ (®i,j+l ■ ®ij )®j 

j-1 


The convergence of radiant flxix in layer i is given by 


■ <+i - ^+1 - < "I 


N+1 

j-i 


j+i ■ ®i+i,j ■ ®i,j+i ®ij)®j • 


(3.11) 


The change in teoperature in layer 1 resulting from the convergence of 


radiant flux is 
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H+1 

S 

J-1 


where 


APjC 

i P 


(®l+l,j-KL ■ *1+1, d 


*l,j+l + *lj) • 


where Ap^ la the pressure difference across layer 1. 

A more complex integration scheme results If we assume that temperature 
varies linearly In log, pressure from layer midpoint to layer midpoint: 


e«) - + a - Y^)e^.;^ 


♦i-i - ♦ - ^ 


where * (An ip ~ Zn 4>^_2^)/(An - An • Above the atmosphere is 

A 

assumed to be isothermal at 0, , and below 4i„ the temperature is extrapolated 

X w 

/N 

from level rf)„ . to a temperature discontinuity at the surface. The expression 
^N“X 

for the Curtis matrix elements becomes 


\j " (*i+l,j+l ®i+l,j " ®i.j+l ®id ) , 


where e^j - e^j + c^j and correction term c^j is given by 


'Id 


- 4 


0, if d - 1 


e , , “ e , , 

id Ij 


'id 


i,d-l 


Yjde, 2 5 j < N-1 


*i,N+l " ®1,N 


"i,N+l 


Yj^de, j N 


1,N-1 

(0, ifd -N+lorN+2 


(3.13) 


The integrals ln(3.L3)have been performed numerically using the trapezoidal 


rule. 


3 , 2 The overlap of spectral llnc« 


Below about 30 km pressure-’broadened line widths may be sufficiently 
large to cause significant overlapping of individual rotational lines. The 
entire band (550-800 cm“‘) Is divided Into ten equal spectral subintervals, 
and spectral lines are assumed to be randomly distributed In position within 
each subinterval. Eqn. (3.5) can be written 





1 10 

H "'i 


0 


k-1 


(3.14) 


where T, ■ 


k 6v, 


dvT is the mean transmission function over subinterval k. 


6v, 


The assumption of random distribution implies (Goody, 1964, p. 154) 




(3.15) 


where W- Is the equivalent line width of the line in subinterval k . 

Xi 


3 . 3 Line shapes 

Through much of the region of interest the shape of the Individual 
spectral lines is given by the Voigt profile, a convolution of the simpler 
Doppler and Lorentz shapes. The equivalent width of a Voigt line can be 
adequately approximated with a formula that interpolates between the 
equivalent line widths given by the Doppler and Lorentz profiles (Rodgers 
and Williams, 1974) 


W(p^,Pj,p) - (W2 + W2 - Wj W|/s2u2)l/2 


(3.16) 
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where 

and 

- Wjj[u/u,S(T^j)l . 


is the Lorentz equivalent line width given by an accurate approximation 

due to Goldman (1968). is the Doppler equivalent line width given by a 

series approximation or a large Su approximation (Goody, 1964, p. 136). 

is the line strength evaluated at the Curtis»Godson temperature 

between pressure levels and and u is the amount of absorber 

between levels p, and p,. The Lorentz line width (p, .>T ) is given by 

i i L ij Ij 


where the Curtis-Godson parameters are 


i/2 




T(p) dp and p^^ =” j (p^ ' 


For a given temperature profile the Cur tis‘~Godson temperature was computed 
for every pair of levels and stored for future referenct^ as a two-dimensiov;..al 
array. 


3.4 Variation of line strength with temperature 

There are nearly 17,000 CO 2 lines tabulated between 350 and 800 cm*“^ 
on the AFGL Atmospheric Absorption Line Parameter Compilation ("McClatchey 
tape**). Rather than sum over all these lines in (3.15), the lines in each 
subinterval were divided into groups according to line strengths, and a 
mean line strength was used to represent the strength of all the lines in 
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that group. Working with one subinterval at a time, the strength of each 
line was computed at seven temperatures between 150 K and 296 K with the 
formula 


S(T) 


S<T5>Q„(T3)Q^(T3) 


Q,(T)q^(T) 


exp 


1.439E” (T-Tg) I 



(3.17) 


where and <4re the vibrational and rotational partition functions 
interpolated from the table in McClatchey (1973). Tg is a reference 

temperature, and E'^ is the energy of the lower state of the vibrational- 
rotational tran.sitioa expressed in S is expressed in the units cm*“^ 

(cm STP)*”^. At each temperature the lines were sorted Into 50 groups 
representing equal increments of logj^QS. The resulting histograms of line 
strengths for one subinterval are shown in Fig. 2. The histograms were 
represented approximately by dividing the lines into two categories: a set 

of temperature-independent strong lines, and a larger set of weaker lines 
which follow the distribution function 


n^(x) - m (x - foe x < x^^ , (3.18) 

where x - log^oS. Here x^^ » where Is the strength corresponding 

to the empty bin of lines just stronger than the strongest temperature 
dependent lines in sub interval i. Factor m wai determined for each tempera- 
ture in each subinterval, but was found to be adequately represented by 
a universal function 


m(T) - aT d- b, a - 0.001615, b 


-0.2066 , 


(3.19) 
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where a and b avc evaluatred by regression. The mean strength in each gt’oup 
is given by 

Sj ■ exp (2. 303 x^) where - Jd:^ti(K)x/Jdxn(x) » (3,20a) 

and the number of lines having this strength is 

Fj J dxn(x) , (3.20b) 

The number and mean strength of the strong lines in each subinterval was 
evaluated directly. These strong lines, together with weak lines evaluated 
from (3.20) were used to compute mean transmission functions, flux 
equivalent widths, and finally, the Curtis matrix elements for a standard 
temperature profile. 


^•5 Non~LTE source function (this section follows Houghton, 1977, p. 58ff) 
Above about 65 km the frequency of collisions is insufficient to 
maintain a population of excited vibrational states consistent with the 
local kinetic temperature. The source function for a particular transition 
is no longer identictil to the Planck function, but is given by 


J 

V 


L -h eB 

V 

1 *f e 


f 


(3,21) 


where e - C/A is the ratio of the rate of collisional deactivation to the 
rate of spontaneous radiative deactivation of the excited state. Assuming 
and B^are constant over the frequency band, and using the equation of 
radiative transfer to eliminate the mean radiance gives 


J • B + q Cp/47TceSj^ , 


(3.22) 


u 


where c is the concentration of the absorber, and and q are the local 
band strength and radiative hei.^ ig of the particular transition. Eqn. (3.12) 
muse be modified to 

N+1 

"l ■ 2 • <3-23) 

j»l 

Using matrix notation, (3.22) and (3,23) can be combined into 


q 


j^nLTE 


0 where R 


nLTE 


(.. - RE)"^R 


and 


E 


ij 


( Cp/4irBce^Sj^^^ ] 


(3.24) 


^ and are the band strength and LTE factor evaluated in the i*"^ layer. 

In the present calculation all lines are considered part of the funda- 
mental band with Einstein coefficient A * 0.74s“^ and collision frequency 
C ^ 0.67 X 10^ s“^ at standard pressure and an assumed mesopause temperature 

of 180K. The band strength S is evaluated as 

b,l 

10 

^b,i 

k-i i 

where and are given by (3.20) and F® and refer to the strong 
lines in subinterval k. 


3.6 Temperature dependence of transmission function 

Since the flux equivalent widths depend on the temperature of the 
intervening material , the Curtis matrix elements are functions of the 


w 
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atmospheric temperature profile. The flux equivalent width is related to 
the transmission function, so we have adapted the Interpolation method of 
Pels and Schwarzlcopf (1981). 

Let k^(p”) - ^f jj^[v,T(p' ') ,p' ' ) where f ^ is the absorption coeffi- 
cient of the line In 6v. Let T(p") - T^(p*') + <5T(r'') where 
is the standard temperature profile. To first order in 6T 


T^(P,P' »W) * exp 






cdp ’ * 

/ ^ 3T 

0 1 vp ; 
T 

gy 


exp 


P a 


cdp ' 


gP 




Assume 6T(p'') is small and expand the second exponent 


T^(p,p' .y) 


T«(p 

V 


L •'p £ 


^ 6T(p")-^^ 
VP ^ 


where 


T^PfP' *y) ” ^*P 


p I 


Eqn. (3,5) becomes 


i(p.P') - J dp u I [l - tJ(p,p',p)| 

® Av 

■/.'■‘ Ft '- L - -.'...•.■>1 


6T(p") 
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i 

Now define «*^(p,p') • J dp p T^(p»P*»^^)j '**'•* 

0 iv 

1 f) f 

M£^,Eia- _|.jdp jdv T0(p,p’,p)y-^ [v,T(p-),p"] ^ , 


t(p,p’) - e°(p,p') + f dp” 

-'n 


5T(p”) 


Consider the functional derivative 6e/6T. Suppose it can be factored 
(and In several approximations It can) as 


"S e(p,p',p”: 

6T 


F(p,p')G(p”), 


Kp.p') - J dp’ ' 


6T(p”) = r dp'*F(p,p’)G(p'')6T(p”) 

J n 


= F(p,p’) f dp”G(p”)6T(p”) 


If we let e (p,p*) represent the flux equivalent width for a temperature 
profile everywhere warmer than by a constant amount A, and similarly 
for e (p,p’)f then from a Taylor series we have for any constant temperature 
deviation 6T 


-I. g-jPt P*) - e 


e(p,p’) = e°(p,p') + 


7|e'^(P.P') - 2e°(p,p') + e~(p 


,P^ 


(3.26) 




17 


Now define the weighted temperature deviation 

P * P ' 

A(p, P’) - f dp''G(p”)6T(p")/f dp*'G(p") . (3.27) 

■^P -^P 

We will assume that for variable temperature deviations (3.26) holds 
with A(p,p’) replacing 6T. The best choice for weighting function G(p) 
and for A have been determined by Pels and Schwarzkopf as 

G(p) « pO-2 [1 + (p/p^)°-®] (3.28) 

where p » 30 mb, and A - 25K. 
o 

Therefore the procedure adopted has been to compute a matrix of e, 

and for each of three temperature profiles* Each R is then 

nLTE 

"unraveled" using (3.12) to obtain a matrix e . This can be accomplished 

because there are "boundary conditions" e^^^ = 0 and e^^^ Av/2. Beginning 

with the lowest level, successive layers of e can be computed from (3.12), 

first on the lower side of the diagonal (j > i), and then on the upper (j i) . 

These three arrays are stored as data for the dynamical model. When the 

nLTE 

radiative algorithm is called, e (T) is ccraputed by (3*26) for the appro- 

nLTE 

priate te-flsperaturc profile, then R (T) is computed from (3.12)* 

The Curtis matrix for the lOy band of CO 2 is computed using the same 
procedure. The only difference is that the lOp band was not divided into 
sublntcrvals to compute the mean transmission functions of (3.14). 


rvRK'- V 
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A. INFRARED COOLING BY OZONE 

Th« version of the algorithm currently incorporated in the dynamical 
model uses s cool-to-space plua cool-to-ground approximation utilizing the 
band absorption formula of Ramanathan (1976). In this section we shall 
develop a more general method from which the above approximation can be 
easily derived. 

Several empirical and semiempirical formulae (e.g., Walshaw, 1957, 
Aida, 1975), have been developed for the integrated absorption A of the 
9.6p band of ozone where 


A 



J dv(l - T^) 
AV 


(4.1) 


These formulae may be used to compute flux equivalent width and Curtis 
matrix elements fro® (3.5) and (3.12). The procedure outlined in section 3 
may be greatly simplified. IR cooling (or heating) by ozone is significant 
only below about 60 km, so non-LTE effects can be ignored. The contribution 
of the 9.6y. band to the total cooling is small enough that variations of 
the transmission function with temperature profile can be neglected. However 
since the vertical distribution of ozone varies, in general, with latitude 
and time, a new Curtis matrix must be computed for every ozone profile. 

Use of the Curtis matrix method facilitates the identification of the 
various components of the radiative interactions. The radiative cooling is 


N+1 


’i ' I! 
1-1 


given by 
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where Hfl !• the ground. If the atmosphere if isothetmal, and the ground 
is at the same temperature^ then radiation exchanged with the atmosphere 
and gtousad provides no net cooling and there is only cooling to space. So 


GTS 

*^1 


N+1 

I 

J-1 




N-hl 





(4.2) 


If the ground is at a different temperature, there can be cooling to both space 
and ground. Then 


CTS 

‘^1 



N+1 


J-1 






^ < Vi 


- 


(4.3) 


vhere C* - 

A further approximation can te Introduced for 0. The ''temperature" is 
given by 

e“ - 1 

9 - ■ — where u - hv/k.T - 4.988 for T “ 300K. 

^uTq/T _ ^ o o 

Hence, e'* - 1 r e'* and for T < T we have e - 1 ~ 


and 


9 = 


so qfS + <1™ = .“c“ + ."c*(.-"^o/Tg . .-oT^/Tij 


. 146.6 Cj + 146.6 - e-^''W/Tl) . (4.4) 
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where T Is given In Kelvins. The radiative cooling of the 9.6vi band for 
the U.S. Standard ozone and temperature profiles computed with the Curtis 
matrix method employing Alda's parameterization Is compared with the cool- 
to-space plus cool-to-ground approximation In Fig. 3. 
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s, progkams and subroutines 

5 * 1 Routines used to obtain line strength distributions 

The 16899 lines attributed to carbon dioxide between 550 and 800 cm”"' 
have been sorted into eight sets according to the energy of the lower 
state of the transition. Each set is written on an element file of the 
flleset C02SET. The description of each element set is as follows ("normal" 
refers to 


File 

name 

Lower 

level 

Number 
of lines 

Type 
of CO2 

FU626 

0 

153 

"normal" 

H0626 

^2 

546 

"normal" 

HT626 

v^*2v2 

159 

n i f 1 

norma X 

HD626 

2v2 

691 


HF626 

3^2 

564 

‘^normal** 

HC626 

i4v2 

3805 

**normal'' 

FUISO 

0 

1269 

.Isotopic 

HTISO 

> 

A| 

9712 

isotopic 


Program BINS sorts all of these lines into bins according to strength 
without regaro to band of origin. The entire band is divided into ten 
equal subintervals and the strength of each line is computed at seven 
different temperatures. Tlie width of each bin is 0.25 in units of log S, 
where S Is expressed in (cm~') (cm STP)"'. An example of the resulting 
histograms is shown in Fig. 2. If the McClatchey tape is used as the input 
for program BINS, then DO loop 12 should be replaced with IB =» 8, the two 
IF statements after statement number 98 should be activated, and LMAX must 
be set equal to or greater than 16899. 
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Frograa FIT fit# each of the seventy histograma ptoduced by BINS 
with a second order ^•'olynomial, according to (3,18). The only free 
parameter la the coefficient m, 

Graphs of m va , T for each of the ten subintervals are shown in 
Fig, A.- It is sufficiently accurate to assume a universal function 
m(T) • aT + b for all subintervals, Program REGRESS provides the least 
square values for a and b. 

Spectral lines associated with the fundamental have no vibrational 
temperature dependence and are usually strong. They can be identified on 
Che histograms because they rise above the generally parabolic shape of the 
line strength distribution profile (see Fig. 2). Program RESID subtracts 
the fitted line strength distribution from the actual distribution, and 
computes the number of strong lines and their mean strength in each bin, and 
then consolidates each four bins to give the number of strong lines and their 
mean strength in each strength decade. From this information one character- 
istic strong line number and strength is chosen for each subinterval, and 
specified in program CURTEL in arrays FS and SS. 

5 . 2 Routines used to obtain Curtis matrices 

Program CURTEL computes the Curtis matrix elements for the 15p band 
of carbon dioxide. Wave numbers are expressed in units of the Doppler 
width at 300K. 

The computation of the non-LTE Curtis matrix begins with the evaluation 
of the total band strength by summing all the mean strengths weighted by the 
number of lines at that strength over all strength decades and subintervals. 
VMULFM and LINVIF are IMSL routines for matrix multiplication and inversion. 
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Subroutine FbUXEL calauXates the equivalent band widths, First the 
pcessute levels (PHI) and layer midpoints (PHIH) are computed. The Curtis- 
Godson temperature (TCG) between any two levels is computed and stored. 

The correction terms c^^ (see (3,13)] »t <i evaluated by the trapezoidal 
rule with the interval divided into Ml oanels for the diagonal terms > 

M panels for the nondiagonal terms, and 3*M/2 panels for the terms at the 
top of the atmosphere. 

In subroutine EFUNC there is an opportunity to use the integrated absorp- 
tion parameterizations of Pollack jMt (1981) or Rodgers and Walshaw (1966) 
instead of the line strength distributions that have been derived. If the 
integrated absorption parameterizations are not chosen, the Curtis-Godson 
temperature and pressure are evaluated and the line strength distribution 
parameter [m in (3.18), called A2 in this subroutine] is calculated, 

Subroutine W performs the summation in (3.14), evaluating the mean 
transmission function for the entire band, The line strength distribution 
parameter is used to compute the number of lines and mean strengths in 
nine decades of line strength for each of the ten subintervals. This 
synthetic line distribution is cut off at the low strength end so that 
the total number of lines in each subinterval is the same for the synthetic 
and real distributions. 

After the set of non-LTE Curtis matrices have been produced for the 

three temperature profiles and for both lOp and 15p bands, program EQWD 

inverts (3.12) to obtain a set of equivalent band widths. The set of six 
nlaXE 

arrays of e are used by subroutine CURT to compute the temperature- 
dependent Curtis matrix elements. 


There are, unfortunately, two different conventions for the labeling 


of the layers of the models In the dynamical model (JRH) the lowest 
stratospheric layer is denoted J ■1, and the highest layer is J ■ JML 
(i.e., 16). In the radiation algorithm (CBL) the highest layek is 
denoted J’ « 1, the lowest stratospheric layer is J * • JML, the tropo- 
sphere is j’ » JML ^ 1, and the ground is j ’ JML 2 (i.e*, 18), 

Comment cards in the listing noting the shift of conventions will hopefull 
prevent confusion. 

^•3 Routines used for ozone and oxygen 

Program OZPRO produces the arrays needed to compute the radiative 
heating and cooling due to ozone, given the ozone mixing ratio and 
temperature profiles. Subroutine COLOZ computes the specific density and 
column abundance of ozone. Subroutine CMOZ96 computes a Curtis matrix 
for the 9.6p band, and from it cool-to-space and cool-to-ground coeffi- 
cients are derived. Two different parameterize tions (Aida, 1975 and 
Walshaw, J957) of the integrated absorption of the 9.6p band are offered 
in function OZ96A. 

Subroutine OXY computes the density and column abundance of molecular 
oxygecc. 

Function 02UV computes the solar energy absorbed per molecule of O 2 
from the parameterization by Strobel (1978). 

Function OZUV computes the derivative of the fraction of total solar 
flux absorbed with respect to ozone column abundance from the parameterize 
tion of Lads and Hansen (1974). 
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5 • 4 Routlnea uaed In th« dynamicaX model 

Subcoucln« HEAT computes the dlabatlc heating term for the dynamical 
equations. The earth-sun geometry factors and the set of Curtis matrices 
are computed once per day. The static stability is computed once, during 
the Initi. 1 time step of the run. The zonal mean heating is Introduced 
into the dynamical equations gradualTy with a growth time (GlIME) of two 
days ♦ 

Subroutine RADEQU computes the global radiative equilibrium temperature 
profile (i.c., the temperature profile such that the radiative heating 
integrated over the entire globe is zero at each level). 

Function DELT computes the individual terms of the net radiative heating. 
The parameterization of Lacis and Hansen (197^i) is uaed to compute the 
radiative heating due to absorption of solar visible and UV radiation by 
ozone. Values of the various parameters describing the shortwave planetary 
albedo are taken from the program of Ramanathan (private communication). 

The parameterization of Strobel (1978) is used for the absorption of solar 
radiation by molecular oxygen, 

5 . 5 Changes in main program WAVE 2 

The flux equivalent width arrays are read in first for lOp, then 15p, 
in the order minus, zero, plus, KAP is the Newtonian cooling coefficient 
(day~^) derived from the present algorithm, assuming a U,S, Standard tempera- 
ture profile and the ozone profile (based on Krueger and Minzner, 1976) 
currently used in the dynamical model. The coefficients have beer, computed 
for a "boxcar" perturbation of vertical scale 25 km. TZO is on input the 
U.S. Standarc!i aj^mospherc temperature profile. Array PRS contains the Fourier 
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(ifeooniposiclon of the observed u at the lower boundary. QOG and QO are 
Q 3 (ground) and Q 3 of (2.1). QOZSG and QOZS are n 3 (ground) and 
of (2,2), CTS la -146, 6C® and GG is 146, 6C* , defined in (4.4). 

Subroutine HEAT la called every IRCT time steps; IRAD Is used as 
a counter. 

Tropospheric and ground temperatures are computed with simple equations 
that fit the climatological temperatures to within a few degrees. 
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Figure 1: Vertical structure of atmospheric model. 
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Figure 3: Comparison of cool--to-space plus cool -to-ground approximation 

(broken line) with Curtis matrix method (solid line) for 
radiative cooling by the 9.6y band of ozone. Temperature and 
ozone profiles are from the 1976 U,S. Standard Atmosphere. 
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Figure A: Weak line distribution parameter m ^ temperature. Solid lines 

connect values for the same subinterval. Identified on right 
(e.g., subinterval 1 is 550 ^ V < 575 cm“^). Broken line is 
best fit for all sublntervals . 









APPENDIX 


FORTRAN Code 


PRUORAH • IHS (INPUT, OOTPUr/TRPei»TAP6<#!Ai?6J*TAi»6^»IAf 15 
*,rAPvt>»TA^£7»TAPk8,r*PSV) 

C SUbROUUNE TO iUhT UlfkS ACCOROINC TO ifSieNCTrt 
C (TAPCS 1 THROUGH 8 *Rt fcltMFNT PUtS OP C02SITI 
C 

OlrtSNSTOH N(10,7,50T 
.OINFNSION TP(7>,OV(7> 

IMTIS(R US 

UiTA TP/ 180., i7>.»200 ,ZSO. ,2 75»,256,/ 

0 4T* OV /l.Oo A, X. 00V&, i.0l9Z, 1.0327, l.v502, 1.0715, 1 ,0</31/ 
TS*Z56. t 0VS*l.O531 \ PJ'l, 

OSTP»2.to0696»19 
C Jflin OUTPUT APRAVS 
00 U IS -I, 10 
DO 5 IT»1,7 

00 J 1*1, pU 
H( IS, IT, n«o 

3 CONTIhUfi 
9 CONTlNUe 
U CONTtNUe 
C SORT EACH BAND 
UO 12 

1 S* I 

UPLI.i«575. 

P£AUUB,15» LNAX 
LL«0 

99 LL*U»l 

IF(LL .GT. LNAX) CO TO 10* 

98 RfcA0(lU,5) FRcU,S,Al.PhA,b(,S,US,LS,ROT,0ATe,ISO,t1OLC 
C iHEHtO .LT. 550. » GO TO 98 

e IF(N0LC .NE. 2> GO TO 99 

IP(FREQ .CT. UPLIH) CO TO lOl 
97 COtiTliSUE 
DU 7 ir-1,7 

^•rp^m 

tF*EX> U.<.39*ELSHT-TS>/(TATS» > 

ST-S»EF*QvS/gV( m*n /rs)**(-RJT 
C tXPPFSS STRuNGTh in units (tn ATN)-! CH-1 

ST»Sr*OSTP 

SLT*AL0G10(ST> 

C CHOnSk CDRkECT BIN 
IG»50 

DO 31 I»1,A9 

P*I-pO i RL-1.25>»/A. 

IFISLT .GT. BL> 1G*50-I 
31 CQMINUe 

C count Ll.thS IN EACH BIN AND STORE PAKArtfcTERS 
N(IS,ir, Iu)*N(iS,IT, IC»*l 

7 CONTINUE 
CO TO 99 

101 CONT.NUE 

IFUB .Nc. d) GO TU 13 
DL-UPLIN-25. 

WRII£(9,3P) uL,UPLin 
00 8 IT- I, 7 
RB1T£(V,A5.> TPIITl 
DO A I»1,P0 
NN.N( IS, IT, I > 

IF(NN .Ew. 0> GO TO A 
WR1T£{9,25) I,KN 
A CONTINUE 
lP-51 

Mr1TE(9,29» ip 

8 CONTINUE 

13 UPLlH-UPU.1-29. 

IS-tS-1 

1F(UP11H .CT. bOO.I CO TO 12 
GO TO 97 
12 CONTINUE 

25 TuRNAT(IOX,2I6,E15.S) 

5 .F0K«AT(F 10 . 3 , E10.3,F5 .3,F10.3, U2, 3 X, 1 12 , A 7, A A , A A, 2 X, I U 
IJ F0R.1AT UOX, I5T 

35 FOftnAT(* iPECTPAl INTERVAL*, Fp . 0 ,* TO*,F5.0> 

A5 FOhHATUOX,* Tc NPE R A TORE - F 5 . 0 ) 

END 


fROCRAM FmiNRUT^OUTHIT»Uf(EUT*fCJ> 

C SObHOUTli'il 10 FIT HIST0C,»4MS VITH FAK*3 UlAS 
C tlAPCl IS oUlFUt OF FROGhArt 5INSI 
C 

0IHCN$I0i4 nisd.lcuoi 
data LC/16*«»B» J,0* 

C tc IS •li* NUHrftR OF WEAKEST EMPTY UN IN EACH SUBINTERVAL 
00 4 

REAi>ti»4$) DUN 
DO J IT-l»T 
REAOllfSI T 
P»l. 

TOT-0. 

00 I I-i»M 
Nm«o 

1 CONTINUE 

DO 2 1-USl 
FEAOII^ISI Il^NUn 
IFUl .C0> ?1) GO TO 6 
TN«N(U I 
TOT-TOr*TN 

IFITN .IT. PI GO TO 2 
NC • i I 

P*TN 

2 CONTINUE 

6 IFdJT .EO. 0.1 GO TO 3 
IFINC .CT. 40> NC-AO 
K)>LC(ISI 
NP-NC-lClAl 

Sl-53-0. 

00 7 l-l»NP 

si*sx*fi(icui-i 

S2»S2*< I-U«*4 

7 CONTINUE 
RM-Sl/32 

WRITfc(2»25) ISUCnS>/T,NC» TOTfRM 

3 CONTINUE 
kRITE (2^3:*) 

A CONTINUE 
3 F0RmAT(23X»F5.0» 

25 FORMAT (2 15/ F 10. 1 , 1 5/ F 10 . 1 , 2 FIO , 5/ E 10. 4 1 
35 FORMATI/I 
45 FORMAT! 10A6I 
15 FOPMAmOX/2I6,El5.5» 

END 


PHOCtA.I REGUSSi ( INf Ul ,OUT?UT> TAP Cl I 
C SUeRljUTIHC ru COH<»UTIc •cst values COA. UNIVCRSAL flf 
c irtfu IS dUTPur oc j'*uokA« cin 
c 

RCAl n#H 

s«-sr 2 »snl»$T»o. 

N*0. 

uo 2 I-UIO 
00 3 If* 1^7 
kCAD<W4>) frlf 
miNf 5» ffd 
K*M*l. 
sT.$r»i 
sni.s«i>H*r 

, 5T2»SI2*T*f 

3 CUMflNUf 

AEAQIltlS) OUN 
KEAd(l#lS> OUA 
2 CONIIKUC 

0»N*S I2-S r*5f 

A-<.’<*!/nT“S'1*Sn/0 

9 *<s«*sr 2 -s>iiAsr>/o 

PRINT 25/ A/ 9 

5 PnRrtAT(luX/Fl0.1/l5X>Fl0.5» 

15 FORrlATdOAbl 

25 FUK.nAT(* A»*/£13.5/10X/* e**/E15.9) 

END 


fftOGRAfl HtS10(iN^Ut#OVTf UT» fAPEll 
SUBROOTIMC fO COJIPUTfc HfSIOUAl STRONG U.leS IN 15 HICRON SANO 
HaPCI tA UOtPOT OF PROCRAI BINS) 

OlHtNSION N<IO#7#9l)»*ZnO>#A(7>i.SK(7 >»S$I7I»«<7I 
DATA AZV16 9*# B ♦# i» > 0* » I *# / “T*# li 16# / 

UO A IS*l»10 
AEAd(i#65) OUN 
60 1 IT- 1*7 
ft£A0U*5| T 
00 1 1-1*51 
N(I**ir*I)-0. 

1 CONTINUC 

A(m-0,OOinB<»PT-e.20l7 
60 2 1-1*51 

AEA0a*l5I lI*NU$*iT*m 
IFUt .EQ. 511 GO TO 3 

2 CONTINUC 

3 CONTINUC 
A CONTINUC 

UO 6 IS -1*10 
PRINT >9 

RtSlOUtS FUR STROnCCST OfcCAUC OF LINt STRgUGTHS 
OU lA TT-1*7 
RK-N(IS*1T*L1 
SSI ITI-0 
R( iri-0. 

IFIl .LC* XZIISI) GO TO 16 

K(iri-RR-AiiT)-n.-xzas) >-*2 

16 SRUTl-ilim 
1—0.25*11.-5.51 
S-10**C 

■IFIkR .CO. 0.1 S-0. 

SSI m-s 
lA continue 

PRINT 35* ISRIITl*IT-i*7} 

PRINT A5, ISSim*IT-l*71 
RCllUUcS FOR SEVEN HORc UcCADES OF LINE STRENGTHS 

00 9 K-l*7 
00 12 J-1,7 
SRI Jl-S$tJl«0. 

12 CONTINUC 
00 d J»1*A 
I-A*IN-1)*J*1 
RI-I 

00 II IT-1,7 
RR-NI IS, IT, I » 

RlITl-0. 

IFIRI ,LE. XZIIS) 1 GO TO 13 
RIin-RR-AIin*IRl-XZIIS »)-*2 

13 SRIin-SRIITHPIIT) 

E--0.25*IRI-5.51 

S-IOR-E 

SS(1T1-SSI1T1*$*RI1T1 
11 CONTINUE 

PRINT 25* I, IRIIIl, 11-1,7) 
e CONTINUE 
00 17 IT-1,7 

IFISRIIT) .EO. 0«1 GO TO 17 
SSIITI-SSIID/SRIIT) 

17 CONTINUE 

rtCAN FOR FOUR BINS lONE DECADE) 

PRINT 3>, ISRIII),IT*1,7) 

PRINT A5, ISSI11),IT-1,7) 

9 CONTINUE 
6 CONTINUE 

25 F0R«AIU5,7F10.1) 

35 F0RNATI7Fl:,.l) 

A3 F0RNAI17E15.3) 

IS F0RNATII0X,ZI6,E15.5) 

35 FORMATIlHl) 

5 FORMAT I23X,F5. 01 
65 FORHAfllOAb) 

ENO 


fHOGRAH CURTEUIM^UT^OUTfUT»TAP£fc»TAft7> 

C0MH0N//S(10.> „0),F( l0a0)»0£lNU» ALPI40Z* El 141*40 )^E(A1^40> 

** PHU40>*PHIH(40)*rCG(20»20)f TT4 2U 
OiMEi^SIQN T(21)»4<41/A0>,0420)*CTS420»*THETA4Z0» 

OiricHSlON H(2I*Z0>fU(20*2O)?UI(20*2O>*UK(2O)*K4U0*21)*QO(ZO) 
OIHENSION FSUO)*S$nc) 

REAL NU2»KB*RASS 

DATA N*rtf ni/L*0ELTAHU/17*4f 16* }*2»0./ 
data F$/2*0.*20.*).*22.*7.*4A0.F 
DATA SS/2*0. *0.05*2.* 5. 5* 4. 2* 4*0. / 

DATA T/ia 7,1* lb 6. 9* 103.2*201.6*210.4*221.0*234.7*249.4*263.5* 
A< 70, 7* 2 5 A. 0*2 43. 2* 22 9.6* 224. 2* 2 10.6* 216. 7* 231 ,4* zee. 2* 3*0. / 
Pl-J. 1415026 
SPl-SOKI (?n 
HP-6.6237E-27 
AVG«6.02Z5E*23 
fUSS»44./AVG 
C-2.9979E*10 
KB-1.3605E-16 
CP»l.E+07 
GRAV9B1. 

T2-300. 

PAIUnETSRS FOR 15 niCRON BAND 
NOZ*b6/.379 

0-EXF (H?»C*NU27 (K5*T2 ))-l. 

B342.2.»HF*«C*NUU**3/(C*0> 

C0NiT-Z."Pl*B3AP*GRAV/CP 

ALPHOZ*NU7*SORT(2.*KB*TZ/nA5S)/C 

DbLNU>OELTANU 

M»N*1 
DO 2 1-1*10 

F<l*l>-FS(I) i SiI*l)-SS il)/ALPH02 

DO 27 J-2*10 

F(i*j»>sn*j)-o. 

27 C06UN0E 
2 CCriTiNUE 
DU 26 I-1*N1 

mn-rm 

26 CUNTINUE 

PRINT 101* N,H*«l*U*OEtNU 
CAbL FLUXEUN*h*Ml*L) 

DO 11 1-1*10 
00 14 J-l*10 
S(l*j>-sn* J)*ALPHOZ 
14 CONTINUE 

WRITE (6* 103 )< (F ( I* Jl*5( 1*0 ))*<>• 1*10) 

13 CONTINUE 

COMPUTE CTE CURTIS MATRIX 
UO 20 I-1*N 
CTSm-O. 

0P-l.Ol3E»06*IPHI(lAl )-PHim ) 

OU 21 J-1*N1 

R( J, I)- tl{ J*l, I*l)+El(J, n-El( JA1*I)-£1U*I*1) 

M U* n- (CONST /OP )ARU*n *66400. 
cTS(n-CTsm*R(j*n 
21 CONTINUE 

20 continue 

WRlTE<b*100)( (R(i*I ), J-l*Nl)*I«l,N) 

WRIT£(t>* lOOMT(I), I-1*N1) 

VR1T£(6*102) 

WRITE(6*100) (CTS(D* I-1*N) 

WR1TE(6*102» 

COMPUTE TENPtRATURE THETA • B/BIREF) 
v:-HP*c*NU2/(Ke*Tn 
00 16 J-1*N1 
V-HP*C*NUI/ (K9*T(J ) ) 

TH£T4(J)-(EXP(VZ)-1. )/(6XP(V)-l. ) 

16 CONTINUE 


r» r% 


C COKPUTi HEATIH6 HATfS 
UO 10 
OtU»C. 

ou 0 

o(n-o(n«Rij^n*THir*(J> 

•P COMUNUE 

UKlTe(6«10»> Q(II 
10 CONYINUE 

COnPUTC NON-LTE CUM IS HATMX 
CN>0.0<03)*4A./20.V66 
SOANO-0. 

C COnPUU rOTAl band sibencth from oisjauutioh 
DO 21 l^XplO 
00 24 J-lflO 

S6AN0>SBANDfF(lr<MPS( I> J| 

24 CONTINUE 
21 CONTINUE 

SB AND- SB AND* 224 14./ 44. 

C 

C PENFUftN NATRiX ALGEBRA 
00 6 J-1#N1 
00 7 K-V,N 

C EPS - COLLISION R A f E/ SP bNTANE OUS ENISSION RATE 
EPS«0.74*PHIH(KI/1.5E-05 

IF(J .EO, K) HU^K)-CP/(4.*Pl*CH*SBAND*EPS*eBAR*66400. ) 

7 CONTINUE 
6 CONTINUE 

CALL VhULPHIP»H,Ni,N>N»41#21>U,20> IER» 

PRINT 101, lER 
00 8 1-1,N 
00 LI J>1,N 
U(I, J)--U<1,J) 

1F(1 .60. J) U(1,J)*UU,J )*1. 

11 continue 

8 CONTINUE 

CALL L1NV1F<U,N,20,UI,0,WK, lER) 

PRINT 101, lER 
00 17 I-1,N 
00 18 K>-l,Nl 
SUH-0. 

00 19 J»1,N 

SUN.SUM*UH1, J)*RIK, J ) 

19 CONTINUE 
RR(K,l)*SUn 

16 CONTINUE 

17 CONTINUE 
RRITEIN, 102) 

C NON LTE CURTIS MATRIX 

VR1TE(6,100M (KHIJ, I) , J«1 ,N1), I*1,N) 

WRITE(7, lOOM (RRIJ, 1), J«1,N1), I*1,N) 

C COMPUTE HEATING RATE FOR NON-LTE CURTIS MATRIX 
00 1 1*1, N 
CQUT-O. 

DO 4 J-1,M1 

COT I )*QO(I)*RR( J, 1 )*THETA( J ) 

4 CONTINUE 

WR1TE(6#105) OOm 
3 CONTINUE 
DO 12 1*1, N 
CTSUI-0. 

DO IS J-1,N1 
CTS(1)>CTS(I)ARRU, I) 

13 CONTINUE 
12 CONTINUE 

WR1T6I6,10C) ICTS(I), I*1,N) 

103 FORMAT(ElO.l) 

103 FORMAT! 10E12. 4) 

102 FORMAT!/) 

100 FC'RrtAT!9FB,3/9F8,l) 

101 F0RMAT!4I3,F6.1) 

STOP 

END 


SUDS^OUTIM^ FlUXIUM^H#?nAl 

C SU^^ROUtlHe. FOR CAlCUl^flHG FLUX tOUIVAlCH? 84HO WIOfHS 

CCHiUN.^/S(lO# 10)#F< 10#10l#DFtNU^AlPNOZ»Fi(4l#<iO)^M41»40l 
S#FHnV0)l#FHiHU0l#rCGt20#^0>/f i2U 
DlflcNSlOW C(AlfiO) 

Nl«N^l 

OU« DctHvi/2. 

Ufa .CO. 21 0€L-tt9« 

C 

00 1<I I«lfN 
PHim- PFUHCa#NI 
10 PKIHm^\FHFUNCn#H> 

I# 

PHlMtNi-SOPKO.ll 

C 

C COflPUTF CURnS-GOOXON TCHPFRATURe FOR EACH PAIR OF lEVElS 
UO ZZ l»l#Nl 
00 23 

TCGU^J >-T<lI 

IFU .to. J) GO TO 23 

SUi-O. 

00 26 

SUrt.SUHMnU^TU^l) )^(PMIHil^n-PMIHU)) 

26 CONTINUE 

TCGa#JH0.5P$Uf1/CPHlH(J)“-PHlM<lll 
TLOU#n-TCG<I# J) 

23 CONTiNUE 

PKInT 101#<TCGU#J)# J-l#Nl| 

22 CONTINUE 
C 

00 70 I-UHl 
Du 50 J«1#N 
IF U.tQ.U GC TO A6 

IF GO TO 30 

IF U^GE.NI GO TO AO 

C 

EFUHC(PHIU)#PHIUI>L} 

COR-O. 

C CONFUTE NONOIAGONAI TERN 
00 15 K»UN 
EKl-K-l 
EK»K 

PHIK* Hh-6Ki)/tt)APHlHI I 

PHIKI* ( (M-EK)/N)*PHIH(>i-lM (eK/rt)^PHIHU I 
15 COR* CORt <GFUNC(PMIK / Fh I H ( J )# PHI H t J-1 )) ♦GFUNC ( PH IK 1 , PfU H U ) 
3,PHIH(J-1) ) )^(EFUNC(PhIKl,PHn i)>l)-EFUNC (PH1K#PHH UrU I 
C 

CU#D- EFUNC(PHIHtJ)^PHI (Hi U-E(J#I)-C0R/2. 

GO TO 50 
C 

C COhPUTH DIAGONAL TERN 
30 CORON-0* 

CORUP-0. 

C 

DO 35 K*l»ni 

EKl-K-l 

tK-K 

PH 1 K*( (Nl-EKl )/«n*PHIH< I-U ^ ( EK WNl > ♦PHI ID 
PMIKl*(^Ml-EK)/rtl l♦PHlH^-UMeK/rtU♦PH^(^ 

CURDN* CORONMGFUNC<PHIK# PH IH (1 ) ^ PH IH iI-1 ) ) ♦GFUNC t PHI K 1 ^ PHI H (1 ) # 
IPHiMi J-IM )^(EFONC(PHlH,PHl(n#U-fcFUHC<PMlKl#PHni )#m 
PHlK*MfU-EKn/fU )*PHU I I ♦< tKi/rt U f PH IH ( II 
PHUl*Uhl-EM/hl)tPHin )♦( EK/rtl )♦PH1H^ I 
35 CORUP* C ORUP^ IGFUNC < PhlK# PH IMU I ^ PH IMil-l) I^GFUNC < PH IK I, 

IPHlHIllfPhlHi I-im*UFUNC<PHlKl»PHmi >U-EFUNC(PHIK>PMUI |#LM 
C 

CU>n* feFUNC<PMlHm/PHim#L>fK0R0N-C0RUP|/2. 

c 

GO TO 50 
C 


C CUMfUTE NOHOUCON/kL TERU AT TOf 
A6 C(l*n>0. 

E(1*II« £FUNC(0.#PHI( n>t.) 

60 TO 50 
C 

At ESH*1I- EfONC(PH!<NI,PHl(n»U 

E(N1*1I> eFUNC(PHl(Nl)»PHl(lt*L) 

COR-0. 
hJ» 3*N/2 
PO A9 
im-K-i 
EK-K 

PMIK- ((«2-EKl»/N2»*PhIH(h-l>»(EKl/M2>*PHHHU 
PHIKA- ( IMZ-EK) /H2» ♦PHI H( K- U HSP / HZ » ♦ PHI ( N1 ) 

A9 COR • COR»IOFUMC<PHIK,PHlHIfO»PHIH(N~I M»CFUNC(PH1KU PHIHINI 
If PHIHIN-n ) )P(F.FUNCIPHIM,PHHU,L)-EFUH*;(PHiK,PHI (I>»L) > 
EINlf I>-E(H.II-COP/Z. 

C 

50 ElUf l)-E(Jf IltCOf I) 

C 

EHMl, I»-EFUNC(PHHNl),PHH 1»»L» 

70 El(NZfIt> DEL 
C 

RETURH 

END 


F UfiCT I OH PFUHCUflU 

« 

FFUKC*O.I*EXP(«JH* j./7i > 
IFU .E3. 1) PFOltC'O. 
kETUAM 
EHO 


FUNCTION PHFUNCUfNI 

C THIS subroutine CONFUTES PREiSURt HALF LEVELS 
PHFUNC-iURT IPFUNCUfN >♦ PFONC U ♦! » N ) » 

IFO ,E0. 11 PHFUNC-O.SPPFUNCIZf Nl 
RETURN 
C ENO 

FUNCTION OFUNC(AfOfC) 

C PLANCK FUNCTION LINEARLY INTERPOLATED IN LOG PRESSURE. 

GFUNC* ALOG(A/C )/AL0Gta/O 
RETURN 
ENO 


FUNCTinK EFUtlC(A*e*L) 

con^toh//s(io«ioi>Fno>io)«crL'iu**iPHoz>fei(^i>^o)*EUi«^oi 

»^FHI(^0WPHIH(4)0t» TCG (20*20 )> (1(21) 

OIMENStON Y(2» 

DATA V/0. 2113*0. 7007/ 

FLEV(XI>7./5.*AL0G(10<,*A)>17. 


IFIL 


u 

CO 

TO 

10 

IFU 

• bQ. 

2) 

GO 

n 

20 

ifii 

• eo. 

1} 

CU 

TO 

30 


C FOUACK 13 HICHON BAND* SL P* lOOOHtt* C •. 00033. 

10 FFbNC«<o.33*AL0GU.»1.13 3b*((9.:6. A3*A6S(A-3n**0.5B6t 
!•( (A3BI**0. 32311 
GO TO 80 
C 

C FOOGEAS-UAUHAW 13 hlCROH BAND* SLP-1000na*C«.U0033 BY VOLUhE 
C EVALUATED AT 200K 

20 EFUNC-0. 

U>ABS(A-B) 

iF(U .EO. 0.) GO TO BO 
P»(A^8>/2. 

TRAN-0. 

DO 21 J-l*2 
UP-U/Y(J> 

TRAN»TRAN*0.5AY (J »*EXP(-300.3*UP/SQRT( l.»5A1.3*UP/PI I 

21 CONTINUE 
EFUNC-170.*(0.5-TRAN> 

GO TO 80 

C 

C 00PPLEP.-L0RENT2 LINE 
30 CRAV'Vei. 

EFUNC-0. 

C COMPUTE PATH LENGTH 

U><B-A|P1.013E«0b*0.00033*22AlA. /<GRAV*26.966 ) 

U-ABS (U) 

IF (U .£0. 0.) GO TO SO 
T2-300. 

C CUNT13>G0uS0N APPROXIMATION 
IA-1 

IF (A ». 1 GO TO 2 

IA-RLEV(A1 

EP>RLEV( AI-IA 

1F(EP .GE. 0.3) IA>1A*1 

2 IB>1 

IMB .EO. 0.) GO TO 3 
le>RLEV«B) 

EP-RLLV(B)-IB 

IF(fcP .Gt. 0.5) IB'In^l 

3 IFUA .GT. IE) IA>18 
IF(IB .GT. 18) lB«ie 
T-TCG(IA*IB) 

C COMPUTE mEAK LINE DISTRIBUTION PARAMETER 
SL>0.(j01613 
BINT--0.2066 
A2»SL*T»P.1NT 
P-(A»8)72. 

P2-1. 

AL-0.08/ALPH0IAP/PZPTZ/T 
C TWO POINT GAUSSIAN QUADRATURE OVtR ANGLES 

TRAN-0. 

00 22 J-1,2 

TRAN-TRAN»0.5*Y(J)*W(U/Y(J»*AL, A2» 

22 CONTINUE 

LFU.HC-0ELNU*(0.3-TRAN) 

C 

BO CONTINUE 
RETURN 
END 




▼ 


FUNCTION 

C this $U6R0UTINk CONFUTES HE*N TRANSNISSlUN FUNCTION 

CONN0N//Sll0^1O)»FaU/lO)#0ELNU/ ALFHD2fEU4W<»0>>E(4)»iV0T 
OlHENSION X2< JOWNUO I 

UATX X2/Ib.#9.^9.>3i/0>>1.»S.>9.^11.»16./ 

OAT* N/b«3* lA07»l676f 2640»3 724»2160»19l<;> 1270. 73 A. *93/ 

C N IS total NUNBER of LINES IN EACH SUBINTERVAL 
C XI IS ONE LESS THAN NUNBER Of FIRST OCCUFIEU BIN 
PI-3.1A1S9Z6 
SFI*S0RT4<>n 
U<0. 

0NU»0tlNU/10. 

C TEN SFECTRAL INTERVALS 
00 o I •If 10 

XN»Xini*T3.AN<n/A2>**(l./3. > 

C, XN N0RNALIZE5 CONTINUOUS DISTRIBUTION FUNCTION TO N 
UI-TXN-5.3>/*. 

U2«(XZn 1-9.5)/*. 

00 8 J’ZbIO 

IFIR .CT. UI) GO TO F, 

IF(R .LT. US) GO TO 8 

VU«*.*R*9.5 

IFIXH .LT. VU) VU-XN 

Vf*.*RA5.S 

IFIXni) .CT. VL) VL-XZ(I) 

C 

C COMPUTc NUMBER Of LINES IN STRENGTH DECADE 

f Tlf J)»A2/3.*((VU-XZT 1))*»3-TVL-Xni) )**3) 

C ESTIMATE MEAN STRENGTH IN DECADE 

XMN-0.25*A2*T TVU-XZU »)♦ ♦ A- i VL-XZ ID )•♦* ) /F U , J ) »X Z ( I ) 
$U» J)«10*A|-0.25*TXNN-5.5) ) 

ST1»J)*S(I»J)/ALPH0Z 
B CONTINUE 
7 CONTINUE 
WU*0. 

C SUM OVER LINE STRENGTH DECADES 
DO * J-If 10 

IFISIlfJ) .EO. 0.) GO TO * 

SU>S(lf J)*U 

C GOLDMAN APPROXIMATION FOR LORENTZ LINE EOUIVALENT WIDTH 
WL-SU/Il.*(0.25*SU/ A)**1.25)**0.* 

IFTSU .GT. *.) GO TO 2 

C SCRIES APPROXlhATIOM FOR DOPPLER LINE EOUIVALENT WIDTH 
WO-WOKSU/SPI ) 

GO TO 3 

C LARGE SU APPROXIMATION FOR DOPPLER LINE EQUIVALENT WIDTH 

2 U2-S0RT I ALOGTSU/SPI ) ) 

W0» 2. •( UZfO. 205 6/ UZ-0. 1 33 S/U2AA3A0. 0070 /U2A AS) 

C INURr-OLATlON FOh VOIGT LINE 

3 wv-sorttwlawlawoawo-wlahlawo*wo/(suasu) ) 
ww-ww»Fn»J)*wv 

* continue 

WI-EXP(-WW*ALPHOZ/ONU) 

WWAWIADNU 
6 CONTINUE 
WW/OELNU 
RETURN 
END 


fUHCnoH UOKWI 

C THIS SOttROUflHC COMPUTIS OOPfL^t UHI iWUlVAlEHT UlOTM FROH SEI^IES 
Pl-3. IA1S^265 
SPl«SQRTtPl I 
NF-l 
SUH«0. 

DO Z N«l,6 

HF«KF^H 

H-NF 

SCH-SUHM-Ult^N/tR^SOpJiUNn 
2 CONTAKUF 
WOl — SPl^SOH 
IFAWOl .Ic. 0«) WOX-1. 

RETURN 

END 


phogran tuwoiiNPvir»uuipi)i,iApes,TAPfci» 

C TriiS iht CJnPurtj tGUlVAL-iil dANO FROH GURUS NATIUCES 

UlHE.HSlUti K< ia»17) f PHiUo>»t:U/# ldUCSU^^ONI2) 
kEAL kB>NU2(2» 

Pl>3.SU:>V26 

HP>b.6237c-27 

C*2.VV7V£klO 

CP»1.£*U7 

6RAy«9di* 

TX-300. 

0HU)-290t A DNtUfZ^a. 

NU2(2>>6fe7.37<» 

NU2U)>979. 

C FOR iO AND 13 niCkOH KANOS 
00 1 N-U2 

i>-eA?(HP*C*NU2(.1»/«KB*T2U-l. 

BbAk« 2 .*HP*tC*riU 2 irt» »**3/iC*0> 

C0iiSr-2.*PI*BBAk*GAAV/CP 

CSr«H»«CUNST 

I COr.UMUE 

N-17 A Nl-tUl A N2>Nf2 A NN-N-l 
C CKL VERTICAL GKIO 
Oii iO 1>1,N 
PHKI )*PFUNC(1/N) 
lU CONTINUE 
P(UtNl»»l. 

C CALCULATc FUNCTION E FOR SIX CUkTIs MATRICtS* 

C K-UU)/ROaOUR«(iUUR-U5)»P.OU:> )/HA(13) 

00 2 U«l»6 

coNsr-c>Tm a DcLnU'Onii) 

iFiU .LT. A) GO TO 3 
LuN$r«CSTi2) A DcLNU>UN(2) 

3 CONTI r«UE 

KEAD(3»25I ( UI I»> J-l/Nl Nl 
OO 29 I>1#N 

OP-l.Ol3EA96*IPHI(I»l»-PHHIU 
00 21 J-l/Nl 

kUU )-0P*iUVUI/(CJNST«cbA9O.> 

21 CONTINUE 
20 CONTINUE 
OU 9 I-UNl 

E (lU l«0. 

9 CONTINUE 
C 

00 A I-l/Nl 
c(N2> U*OELNU/Z. 

A CUN 1 1 HUE 
C 

00 7 JJ*1«N1 

J-N2-JJ 

JM-J-1 

00 K II«1#JN 
1-J-II 

E(J>I)*H(J#n-£(Jfl/l»l)fE(J*^l»I)*^EU>UU 
a CONTINUE 
7 CONTINUE 
00 U KK*1#N 
k 'KK-l 
NK-N-K 

00 12 L-1#NK 

I-Nl-L 

J«I-k 

EU/1 ♦ll»»RUU)U< JU/IUUhU » ll-EU*l»l> 

12 CONfINUc 

II CONTINUE 

Rkiri<6#9)l(E(JflUJ«l,N2Ul«l»Nll 
00 39 I-l/N 

0P«1.013E»06*( PHI ( I Alt -PHI ( n I 
UU it. V-1»N1 

K ( jU )•£< JAl» lAl) A£U, n-EU»l, n-6( 
k UU UKIV# l)«C0NSr«86A00./0P 
31 CONTINUE 
30 CONTINUE 



C I>A1M A TJ CHICK 

KfUMT I'it ((K(J*n>4-Uie»#[-i. i7t 
Z CUNTIHUfe 

;>‘ FUKrt*T(i(&C20.13/l*«£20.13> 

2» 7aHnAT(9f().3t 
k'HO- 


f 





SJ»»OUriN£ CU(^T< rr^KN# TST»k> CUM 

C SUttkuUriKe TO caf1»'urt r£i 1 kb»»ArU.*l- 0 k(*LM)i,HT CURTIS HiTRICeS CURT 

C (T1 iS ACTUAL TbHRi.RA TOHc. RROFUV> TST 1* SlANOARO TcKRCRATURE 
C RkaFTLET 

C CURT 

COknON/CUMTlS/»>.imin>f>iUH< iUT* CUi 7 )«li 2 ( i4>«E(6#19>l3l CURT 

»*CU))|1ArU(>*i7»19T CUKT 

OIMMSION T(Ie»WrsU6), CLTZti'y* lO>*TT(i<i)>TSTU-')»R(2rl0>l7T* CURT 

»U(i0«ia)*u(lb»lO)<aTatil CURT 

KEaL Ktt*HUT(2T CURT 

RNL»Ah-l i 1 A2»<tT»2 CURT 

C CONVERT TO CdL VERTICAL CUOKOlNATE CUNVENTrON • CURT 

00 2 R>l»Kr«L , FERS 

TS(rt) •TST(KN-RI FER1 

T «K)-TTtRk-kT FEdI 

2 CONTtlTUE FEB3 

00 a K-RH»kR FE8J 

TS(iU>liT(R> FEBT 

KRMTTOO FE8J 

3 COkTirTUE FES3 

C COdRUTL TfekFcRATURE OEVIATION FR3H iTANOARl) CURT 

UO 1 CURT 

0T<1» *T(1 »-TSU » CURT 

u(t»i»-OTU.» CURT 

3<l»l)-i. CURT 

I CUNTlNUt CURT 

C CJnRoTc AclCHTtO OtVlAUUtV CURT 

uU 3 CURT 

iP-l»X CURT 

OU A J-IF»KF CURT 

D(I,J)*0TU)*Gi(l)*i>TU»*02Ta) CURT 

dU>J)*GlU)»aZ(V) CURT 

IriJ •t'J» RP) D( i » J ) *0 C i» J ) »U.j * (DT C Kr i -0 T < Kk> ) *Ci tKN ) CURT 

Ll»Iti % L2*J-l CURT 

IFUi .GT. L2> OJ TO 7 CURT 

DO b L*Li/L2 CURT 

D«l/J)*OU/4»*OT(U»«v>ia)AG2<U )> CURT 

eu#j t*3(I>V)^GX(L)iv>2IU CURT 

b COHTINUE CURT 

7 COiUlMUc CURT 

A CONTINUE CURT 

3 CONTTNUt CURT 

c nornalue deviation curt 

00 9 1-i/AN CURT 

IP-I»1 CURT 

00 U J-IPfKP CURT 

D(j>n«oa» j)-o( curt 

II CONTINUE CURT 

9 continue CURT 

C CONPUTt rEIIPERATURE-OEPEHOLNT EOUiVALENT (3AltO> RIDTHS CURT 

00 12 rt-1/2 CURT 

oO U I-l«Kr CURT 

DO lA V1/R2 CURT 

tE(N# J» 1) •c(3*t1-l» J/ !»♦ IE <3*rt#J » I »-t (3vN-2* J»1 ) ) ♦ 0 U » I) / 50. t CURT 

iL..S*(c(3»h#V/I)*L(3*N-2»J#I>-2. •Eti»«-1» J/I))»OU#nP CURT 

SD(J»n/b25. CURT 

lA continue curt 

13 CONTINUE ■ CURT 

U CONTINUE ' CURT 

C COhPUTe CUKTIS MATRICES CURT 

PI-3.1A15R^65A CURT 

MP»b.b2J7k-27 CURT 

C-2.9979E»10 CURT 

K3>1. 30U5E-16 CURT 

CP»I.E»07 CURT 

oRAV*931. curt 

T2-300. CURT 

NU2l2>*ob7.379 CURT 

Numi-930* CURT 

OU 19 11-1^2 CURT 

00*E'XPINP*C*N02(N)/(K3»IZn-l. CURT 

a3AR«2.*MP*IC*NU2<,nj (•*3 /(CaOOT CURT 

CUN5T«2.*PI*ttdAA»0RAV/CP CURT 




Ol) If l-l,Kh 

0P-i.Ol3f»06*(»>tUU>U-PHHU» 
00 id J-r,KP 


j/l r-EC(H« 

M«^J#U»JCUM$ r/OP)*R<M,J> U*»S voo. 
l(t CONfiNUlt 
17 COfiriNUC 
19 COimNUlc 
KCfUMN 
iHO 


CURT 

CURT 

CURT 

CURT 

CURT 

CURT 

CURT 

CURT 

CURT 

CURT 


SUdROUTlNt G(KN) C 

C SUJKUUriHt 10 COMPUTE dcIGHTlNG FoUCTION for rCMPtRATURE- G 

C OtKENUfcUT (fUKTlS iIATkIX CALCUlAHOH ' G 

C G 

COnnON/CURriS/PHHU)»FHlH( V>1U7) /G2<Id)/e(6/19«18) G 

l,K(3t>/17/ld) G 

GC(Pl»P2)-5./6,* (P2*M6./5.)-Pl»*(<>./5.) ) ♦(p2»P2-Pl*PU/0, 12098 C 

ANL»AH-1 t Krt2*/U4-2 G 

C CdL VtKTlCAL COORDUiATE LABEL COIVbttriON G 

Pril(i)«0«. G 

00 1 A«2/KM G 

kJn-R-KN 6 

PrtHA)«0.l*£XP(fcJN*:>./7.» 6 

1 COhTlHUE G 

PriUM4tU*PHIH<K/4»ll-l. G 

PH1H11»*0.5*PH1(2) G 

00 2 A>2>KN G 

Pi4iHtK»*SURnPHUA»*PHl (KflM G 

2 CUNIiNUE G 

UU 3 K«1»M< G 

Gl(K)-OG<PHIrt(K>,PHnK»l>» G 

G2<K)-GG(?H1(R)/PHIH(K n G 

3 CONTINUE G 

G24ANT^1)>0. G 

return G 

END 6 


f«OC*A« QrP*OUMPgT#OOTPyr»TAftl»TA?tll 

»/PHi (^0)»PHIh(^0>,TCGUO/20)»0(l8bOlS<m 
OlNfcHSION 0UCU8>>00< 10* 17) / QO^SC (16 Q02S < 18* I 7 ># CT$ ( 18* 171* 
»cTC(ie*i7l 

tl«EN»lOM b(18>*C$(17)*CO(171 
JH-19 8 KN-17 

» KNL-KN-l 

KNP«KN*l 

C COMfUTt OZONt P8dArtCTfc*8S PD« EACH LATUUOl 70H6 
00 2 J«l*l 

ftEAO(l,)MW(Kl*K«l*KNP) 

PEAPI l*15MI(K»*K-l,KhPI 

C U l)k riAS3> MIXING 8ATI0 IN C/G AT PRESSURE HALF LEVELS 

C NUKaEREO OUUf.WARO FROM TOP OF ATMOSPHcRE. T IS TENPERATURE AT SANE L 

CALL GLOJE(KK*W*I*CS*CC) 
bur,(j ).0(ASR1 
CUIiGtJI-oesIKNPI 
0U(J,KN)-Q(KN) 

OOZS( J*KN)«0ZS(KN1 
CTS( J*KH)-CS(XN> 

CTG( J,KN)*CC(KNi 
OU i K«l*KNl 
00(J*R)«0(KN-(O 
OOZS( J*K)-OZSIKN-KI 
CTS( J*K)*CS(KM-K» 

CTG(J*K)-CC(KN-K) 
i CONTINUE 
Z CONTINUE 

WRITc(2*195MOOC( J)* J-1»JNL) 

WRl TE(2.l95HCO(J,KN)*J-l*JNL» 

VtUTE(£*19^ ) ( (0n( J*K)*J-1*ANL)*K-1>KNL» 

URITEIZ* 195M0UZSG( J)*S-1,JNL) 

WR1TEC2*1V5M00ZSI J*Kh)*J-l,JML» 

•.P*TE(2* 195M (OOZS(J*K), J.1,JMU*K-1,KNL» 

•»ftIT£(2*i95 MCTS (J,KN)* J-1* JNL) 

VR1T£(2*195M (C1S( J,K )*J-1, JNL)*K-1*KNL) 

WRITEI2*195HCTG( J*KN)/ J*l* JNLI 
WP.IT£(2,193H ICTG( J*K>*J.l, JNL)*K»I.,KNU 
195 FORnAI (9E10.3) 

5 F0(MATIE-*.2) 

1,1 F0RNAT(F9.2» 

ENO 


SU8R0OTINE GL08c (KN* V*T,CTS*CTG» 

D1 HENS I UN T(16)*CM| Al*A0),Un.6)*CTS(20)*CTG(20» 
COnNaN//EUAl*AO),E( Al* AO) 

A,PHI ( AO)*PHIH(AO)*TCC(20,20>,OOHO),OOZSI 16) 
KNP«KN*l 
DO 70 I»1»KH 
PHI(*)> PFUNC(1*KN> 

70 CONTINUE 
PHKKNPI- 1. 

DO 71 I-1*KN 
PHlrtm-PHFUNCII*PHll 

71 CONTINUE 
PMIH(KNP)»1. 

CALL COLUZ(KN*W,T) 

00 n L»1*KNP 

PRINT 23* PHIIL),PHIHIU,W(L»/00(U*QOZS(L) 

72 CONTINUE 

CALL CN0Z96(KN,CN*CTS*CTG| 

23 FOKNATISEIS.ST 
5 EnRNAT(E9.2) 

13 E0PNAT(9F«.3» 

RETURN 

ENO 


EVELS 


SUartOUTlHE COLOl(KN»U^ri 

i,PHU^o)>PHlHU0)^TCCU0#20)»Q0M9)»O0lS<l«t 
oiHcnsinN v(ie)>mai 
kEAL H3*HAIK 
KNP>KHf| 

HAIft>28.966 

AVC>A.0223C»t) 

h3*AA./AV« 

kM»8.nA3£»07/«Am 

0STI*>2.6a69E»19 

AT-IOS-l.OiaC^OA 

c>9ai. 

ASSUME OEOMc SCALE HE tCHT IS H ABOVE TOf HALE LEVEL 
h«A. iEtOS 

WOU>-PHIHU>*AIMOS*Wm/ U13*OSTP*RA*mn 
yOZS( U -OOi n*H 
00 2 K«2/KHP 

00(K)-PHlH(M*AI«0S*tatK>/(M3*DSTP*EA*T(K» I 
ASSUME b VARIES LlNcARLY WITH PRESSURE 

CuZSlK)-a0ZS(K-U+0.5/(h3*0STP*C»*<W(K-l>fWTK) )*(PH1H(K)- 
♦pHIH(K-1>)p*TM0S 
2 COMTIHUE 
RETURN 
END 

EUNCTION PEUNCN,NI 

THIS Subroutine computes pressure levels 

RJN*J-H 

PEUNC-0.L*EXP<RJN*5.Z7.» 

1P(J .EO. 1) PFUNC-0. 

return 

ENO 

EUNCTTON PHEUNCU»PHI> 

THIS SU6,<UUTiNc COIIPUTES PRESSURE HALE LEVELS 
U'HcNSiON PHIUBT 
PhEONC-SQRTIPHnH*PHlUM> » 

IHJ .cQ. IT PHPUItC»0.i*PHI (2T 

RETURN 

ENO 


FUNCTION GFUNC<A/B»C> 

PLaMCK FUNCTION LINEARLY INTERPOLATED IN LOG PRESSURE. 
GFUNC- ALOG(A/C I /ALOG (8/C I 
RETURN ' 


SUOi<OUnH? CM0Z96<N»B,CS,CC> 

DiheMsrnh tui >>r< ' ll# ‘.o)^o(zo>#ci$( 2o>#ri«eTAUo)#cs<2o)#cc<2o> 
l#CT5<2v)>,SCUOl 
C0rt!10h//6HA1»'.0)>E<M# AO) 

»,PH1 (AO), ?1I»H)*0), TCC (20# 20 )#OOU8)»002SU9) 
klAL KUZ,K&#KASS 
DATA ft#)U/l#l/ 

i#ATA T7187. I#ld6.9# 193. 2 , 201 . 6# 2 10 . A# 22 1 . 0# 23 A. 7# 2A8 . A# 263. 5# 
t2 70. ?,25t».0# 2 A3. 2# 229, 6# 2 2 A, 2# 2 19. 6, 2 16, 7# 2 31. A# 28 8. 2, 3*0./ 

C WchRBcl#^ 4n£riO)>.EN(S 
PI-3. IA1!>926 
iPl-SOftr(Pl) 

MP-6. 62376-27 
M-N*l 

AVC-6. 02256*23 
hA$S-AA./AVC 
C-2,V'}796*10 
KA-l. 38036-10 
CP-l.e*07 
GHAV-981 , 

rz-joo. 

C PARAMcIEKS POR 9.6 MICRON BAND 
MU2 • lOAO. 

0-6XP (hP*C*NU2/ (K9*TZ)»-l. 

88AK*i, *HP- (C-NUZ ) *»3 /(C*0) 

CCt(ST-2, -PI »B8AR»CRAV/CP 
C 

c 

C COMPUTE CURTIS-GOOSON TEMPERATURE FOR EACH PAIR OF LEVELS 

00 72 1*1, NI 
Ou 73 J-I,N1 
TCG(I#J)-TI 1) 

1 F ( I ,L0. J ) CO TO 73 
SUM-0. 

J-l. J-l 

',0 76 L-I,JM1 

.UM-SUM* { i (L) *T (L*l ) )*(PHiH(L* 1 )“PH1H(1 ) ) 

7o CONIINUE 

ICO( 1, J ) •0.5»SUM/( PHIH( J )-PHlH( I J ) 

(vC( J, I ) -reel 1, J ) 

73 CONTINUE 

fRINI 100#(TCC(I,J),J.1,N1) 

72 CONTINUE 

PRINT 101, M,M1 
CALL FLJXEL(M,M,M1) 

C 

C COMPUTE LTE CURTIS MATRIX 
DO 20 I-l,N 
CSID-O. 

uP-l.0l3t*06*(PHI (Ifl)-PHK 1) ) 

DO 21 J-1#N1 

Ru#i>- ti(j*i,i*i)*Ei(j,n-Ei(j*i,n-6i(j#i*i) 

HU#1)-(C0NST/0P)*R(4#I)»66A00. 
tS( 1)-CS( i)*R(J,I) 

21 CONTINUE 
20 CONTINUE 
C 

PRINT 100,( (R(3,I), J-1,N1 )# I-l#N) 

PRINT 102 

print 100 ,(T(n#I»l#Nl) 

PRINT 102 

PRINT lOU, (CS( 1 ), 1-1#N> 

PRINT 102 

C COMPUTE COOLING RATES 

vz-HP*c*Nuz/(Ka*rz) 

00 16 J-1#NI 
VMP-C-rtUZ/ (K8*I ( J) ) 

IHETA(J)«(EXP(VZ)-1.)/(£XP(V)-1.) 

16 CONTINUE 

PRINT 105, (THETA(J), J-1,N1) 

PRINT 102 
DO il J-1,N 
CI5<4)-CS( J)*THETAU) 

CG(J)«P(hl, J) 

CIGI J )-CG( J)* (THETAIND-THETAI J) ) 


SGU)«CTSiJ)»CTGiJ» 

cs ( n«<o* V. >*cs (I > 

CG(ll«(0*l.t*CGU) 

1 \ CGIiTlNtiC 
00 10 I«l*N 
0 ( ll« 0 . 
no 4 

GHl-Q<U«RU»n»THETA<J) 

9 CONTINUE 

pfUNf ioo» cTsu)*cTcm«SGm,o(n 

10 CONTINUE 
C 99 CONTINUE 
C 106 fOP.IATlfO.lT 
105 fORilAriElO.3) 

103 FnK'UnNl2X*E10t^»2X#El0.5»exn 
lOV FOXiUT(8E6.0» 

110 FnR'lATiaEo.3 ) 

C lOs F0HilATII5» 

102 EQRKXII/I 
C 107 F0RIUTI9E10.31 

100 f llF.nAI I9FU. 3/9E6.3) 

'08 FOhNM (9E ;0. 3/0010.3* 

101 FuRrtATI<il 3 ,F 6 .n 
RE TURN 

ENO 


SU8R0UTINE FLUXEL (N» N>Hl ) 

C SueROUUNE FOR CALCULATING CURTIS CUEFFICIENTS FOR NET UPWARD FLUX 

01 MEN SI ON ClU»AO* 

CCjH10N//EHA1,A0»»EIA1,A0> 

».PHUAOI,PMIH1AO»» ICGI20/ 20»,Q0( 10),QO2S( 18* 

M-N*l 

h2-N*2 

C 

C 

00 70 I-l,Nl 
00 50 J«l,N 
IF U.cO.n GU TO A6 
IF U.EQ. I I GO TO 30 
IF (J.OE.N) on TO A8 
C 

EIJ#I»« cFUNClPHIIJ),PHim > 

CflR»0. 

C 

DO 15 K-l»N 

LKl-K-l 

EK.K 

PH IK- UM-£Kn/N)*PHIHI J-1) ♦! EKl / H)*PH1HI J) 

PH IK I- IM-tKI/NItPHlHIJ-llMcK/ni’PHlHU) 

15 CijR- COP - ( GFUNC IPrt IK / PHIH I J »>PHlHU-l))-GFONCIPHlKl/PHIHiJI 
l» PrtIH|J-l ) ) )♦ ItfUNC IPhlM / PHI ( 1 ) )-cFU,',C ( PHIK./ PH, m I I 
C 

C(J*n- EFUnCIPHIHI JI/PHUI ))-EI J/I>-CQR/2. 

GO TO 51 
C 

30 CORON-0. 

CURUP-0, 

Ell/D-O. 

C 

00 35 K-l/Nl 

EM-K-l 

EK-K 

PHIR-UNl-EKU/fU)*PHIHn-U*ltKl/Nl )-PHI 11* 
PHIKl-lll11-cK)/nl)*PHIHII-l)flEK/hl)-PHIin 

COLON- COftOHMOFUNC (PHIK/PHIHI I ( / PH IN I I - 1 )) *00 LNC 1 PH IKl / PHI H 1 1 » , 
I^MIHI I-l n ) ♦( CFUNL IPHIK/PHI I I I ) -cFUMC t PH I K 1 , PH II I ) II 

FHIK-Hf1l-EKl)/nl)*PHimHeKl/fU)»PHUUI ) 

PHlKl-Hr.l-£K*/Nn*PHI| n»| EK/ni )*PMIHI n 
35 CDRUP- CORUPt (GFUNC IPHIK/PH IH) 1 ) .PHIHtI-1 n -GFUNCIPHIKU 
^ IPHIHI II/PHIHIl imMLFUNCI PHIM/PHII I ) I - 1 F UNC I PH I K , PH 1 1 I) > » 

^ CII/II- EFUMCI?HTHm/PHIimMCQR0N-C0RUP)/2. 

CO TO 51 
C 


4 * ca*ft>o. 

e«i*n» ErONc<o.>PHi(iii 

cu TO 91 


4« C(N«tl- €FUNC(PHI(M*PHI(m 

EFUhC(PHUNU^PHl< m 

COR«(>. 
hi- J*h/i 

00 49 K-l*ni 
tK -!4 

PhtK- n! 1 ?- 6 Kl» /l 1 i)**HIH(N-U* (EKl/K 2 >*PHnNU 
PHIKl- « iMi-eKJ/rtiJP^-himh- n ♦ (fK/niJ^PHUNU 
49 C 0 « • C 0 RMG‘=IINC(PIUK,PHIH(N),PHIH(N“ 1 ) »tGFUNC(PHIKl/PHIH«K) 

1 « > >♦ (fcFUHC < PHIKl *PHU ! I )-tPUNC i PHIK, PHI (I >> > 

c(M*i»- E(Hi, n-giH, n-coK/i. 
c 

91 El(J«n-kU»l)PC(J,l» 

90 cor* r Wide 

c 

ti(Ni»n-eFUKciPHiiNn,PHi( m 

fcl(HZ,n- EFUHC(-l.,-l. » 

70 COhTiNOE 
C 

PRINT lOi 

PRINT lOl, I (6 I n» J-I..N )» EI(N 1 » n» El (N 2 » n. I- 1 ,N 1 » 

101 nJPlAT('<P8.3/10Ftt.3> 

PRINT 102 

100 FORNATI IOEIO. 3 ) 

102 rORHATt/l 
RETURN 
ENO 


FUNCTION EFUNC(A»B) 

C THIS SUEROUTINE CurtPUTES FLUX EQUIVALENT WIDTH 
C Note IriAT MANE OF A6S0RPII0N FUNCTION APPEARS TWICE 
01P.ENSI0M Y(2» 

CGNri.)N//El<41/40»»E<41»40» 
i>PHI(40l>PHIM<40)/TCG(20»20t/QO(LR)>002S(18t 
LATA Y/0.2113,0. /e«7/ 

EFUNC-0, 

IFIA .GE. 0.) GO TO 4 
tFUMC-0Z96Al-l. »0./0, » 
return 
4 CONTINUE 

C lA IS PRESSURE LEVEL CLOSEST TO A, lAU IS LEVEL JUST ABOVE A 
lA-lTlA) 
lAU- IU( Al 

C SinlLARLY FOR IB AND I9U 
Id-I T(S) 

UU-IUIB) 

AU-OQZSI lAUU (A-PHIHllAU) ) / I PH I H 1 1 AU» 1 1 -PH I H II AU )) ♦ ( 00 Z S 1 1 AU ♦ 1 ) 
tCiiZSI lAUn 

BU-QOZil IBU»»(B-PHIHl IBU M / I PH I H I I B U» I ) -PH I H I IBU )l ♦ I OOZ S SI BU A 1 ) 
9L0Z3 I 19UM 
U-ABSIAU-BU) 

IF (U .EQ. 0. I GO TO BO 
C CUPTU-GOOSUU APPROXIMATION 
T-TCGI lA/ IB» 

P-IA*BI/2. 

C OUAOHaTURc over ANGLES 
C TWU POINT GAUSSIAN 

CO 22 L-l>2 

EPUNC-EFUNC»0.9*YIL )*0Z96A(U/Y (L »>T,PI 
22 CONTINUE 
C 

EO CONTINUE 
RETURN 
ENO 


FUNCTION 1U(X) 

C THIS SUU^UUriNF CONPUT6S HALF LEVEL NUKHEII JUST ABOVE PBESSUKE % 
IU»l 

IF(X .LE. 0.) RETURN 

IFIX .UT. 0.1) GO TO t 

IU»7./5.*AL0C(10./1.A292*X)A17. 

IFIIU .LT. 1) IU*1 
RETURN 
2 IU>16 

)F(X .GT. 0.31623) IU-17 

RETURN 

ENO 

FUNCTION ITU) 

C THIS SUBROUTINE CONFUTES HALF LEVEL NUNBER NEAREST PRESSURE X 
C FOR CUKTIS-GODSCN TLMFEhATUHE APPROXIrtATlON 
IT«1 

IF(X .LE. 0.) RETURN 
IFIX .CT. 0.1 ) CO TO 2 
R-7./l>.*ALbGU0./l.A2 92*X)^17. 

IT-R 

E-R-IT 

lF<c .CT. 0.5( IT-IT»1 
IFUT .LT. 1) IT«1 
2 IU-17 

1F(X .GT. 0.5623) IT-IS 

RETURN 

ENO 


FONCTIOH 0196A<U#T»f| 

THIS SUBkOUriNE CQ^PUTES INTEGRaTEO ABSORPTION OF 9,6 HlCRON 
02UNE BAND USiUO FOKflUlA OF AlOA (197S) JOSRT PP. 369*403* 

1F(U .OF* 0*) GO TO 6 
029oA*260*/2* 

RETURN 
6 CONTINUE 

Pl«3.1Ai>92654 

S6-393.9 

CN-X CNAfN-X AT 293.2 K 
A*B»0* 

0B«260* 

X-SB^U/OB 

G-0.07B^P^iT/293.)^^C-0.7> 

0-0. 1 

BI«2.^PMG/0 

y«2.^X/BT 

fe»X/SOFT(l*>YI 

029oA«0«4a.-EXPi-E>l 

Z-Al0010(X/2.62) 
nooiFlCATION TO FORNUIA ST UNM 
Z-A3;<Z) 

F-AIOGIO(P) 

IF<U *07. Z.) GO TO I 

A-0. 16^24 ♦FA ♦ 4tO. 4862 3*F*^3 ♦0*54604 AF^*2-0.11218AF ♦0.04140 

A-IOAAA 

B» (-0.0041 2* F4A3-0.009044F*A2AO. 14160 »F MZ-O. 04420* FAa 3-0 .0 4262P 
3F*A2*0f,50156*F*2. 1600 
GO TO 2 

1 A--0.03097*F**3*0.11180*F**2-0.21613*F*0.21748 

. A-10**A 

B*-0.4800AF*Z^1 .720 

2 C--0.01265*F**4-0.07222*F**3-0.107ia*F**2*0.04269*F^0.55000 

CF*1.-C*EXP(-(Z**B)/AI 

0Z96A-0Z96A*CF 

RETURN 


ENTRY QZ96W 

THIS SUBROUTINE COMPUTES INTEGRATED ABSORPTION OF 9.6 MICRON 
0ZCN6 BAND USiNG FORMULA OF WALSMAW (1957) QjRMS 83# PP. 315-321 

0 

IF(U .G£. 0.) GO TO 7 
Qi9bA»138./2. 

RETURN 
7 CONTINUE 

X-( 1. ♦0.102 5*U )/(!.♦ 1.61*U) 

IF(U .IE. 0.1) GO TO 4 
X«0.964«10A*(-0.53*U) 

IF(U .IE. 0.4) GO TO 4 
X-0.317*U**(-0.74) 

4 A-2.11 

G-U*(XA*A)/(760.*P» 

GL-AIOOIO(G) 

TABULATED FUNCTION FITTED WITH STRAIGHT LINE SEGMENTS 

£• 1 . 

1F(GL .LT* -4.2) GO TO 3 
I —0.01 7 6* (GL A3. 8) 40. 992 
IF(GL .LT. -2.6) GO TO 3 
E»^0.1348*(GL^3.0)^0.94B 
3 F«i.l65/SQRT(l.^734.*G)*€ 

Cl96A-X3e.*(l.-10**(-U*X*F) ) 

RETURN 
END 


0%i 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 

OXY 


FONCTIOH OZbMiin 

C :>vidP.OunHc fJ CALCULATE 5UCAR ttilkO^ AbSUKBcO bY MOteCUlAR OXyOEH 

C RARArtErtKiZAriuil BY STAUbU (19Zi> d olOPHY:> Rtb 63 RR. 622:?*-6Z3») 

C 

C SCHUrtAHN-RUhOE CONTINUUM/ 1250-15 ZO A 
tXl«0* 

fcl»-l*E-17»U 

iE(Ei #lT. -lOO.) GO TO 3 
LXl-tXRiEU 
3 Ol-l.it-17^CXl 

C SCUUrtAliN-AONGt COUTIHUUN/ 1520-1/50 A 
EX2*0# 

t2— 2#9E-l9tU 

lEiE2 .IT* -lOO*) GO TO S 

£X2-EARiE2) 

9 EX3-0. 

E 3 — i* 5 AE-i 8 *U 

lf(c3 *LT. -lOO.) GO TO 6 

fcXi-EXP(£3> 

6 

t 9 — i.ie- 17^0 

lF(b9 .LT. -IOOp) go to 7 

fcXA-cXPU^> 

7 U2-(3.A3^LX2-2.0d*EXi-i.35AEXAI /U 
C SCHUf1ANN-KUMG£ BANDS 
0i«9.03£-19 

IF<U .LT. l*Cn3) GO TO 2 

Qi*i . /(0.193<U^9.6AEA0U*:iURrtUI ) 

2 J2UV^itU2iQi 
RETURN 
END 

FUNCTION OZUVUiT 

C SUBkUUIIHE to calculate SOlAk ENERGY ABSORAgO 8Y OZONE 
C COEKIVITIVE OF LACIS-HANSEN EXPRESSION) 

c 

FI • 1.0^(0*042AU)t(3.23E-A^<U^^2.0)) 

F2 • (0*0212/Fl)bI 1 < U/F !)♦ (OiOA2A <6*A6E-A AU) ) ) ) 

FI • 1.0Am6*6AU) 

F2 • F2t(l*082 / (Fl^ *0.805 ))'*(1.0-U138.6*0«305*U)/fl)) 

FI • 1.04 ni03.6*U) **3.0) 

OZUV * F24((0.0658/Fl)*U^0-<3.0*4Fl-1.0)/Fin) 

RETURN 

END 


02UV 

02UV 

02UV 

02 UV 

J2UV 

Q2UV 

J 2 UY 

02UV 

Q 2 UV 

02 UV 

U2UV 

02UV 

02 UV 

02UV 

02UV 

02UV 

02UV 

02UV 

J2UV 

02UV 

J2UV 

U2UV 

J2UV 

02UV 

J2UV 

02UV 

J2UV 

02 UV 

02 UY 

J2UV 

02UY 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 

OZUV 


SUDROUTINE OXYIT/KNI 

C SUBROUTINE TO CONFUTE DENSITY AND COlUHN ABUNDANCE OF 
C NDIECUIAR OXYGEN 
€ 

C ASSUME N(02)*0.2095*NU0TAU TO 155 KM/ AND N(02)-l.E4)7*EXFI-Z 
C 5/30 KM) ABOVE 155 KM 

C 0 Mrt 0 N/ 02 /D 02 a 7 )/CA 2 a 71 
DIMENSION m?) 

REAL KB 
KNl«KN-l 
Ki*l#380?E-l6 
DO 2 J*1/KN 
RJN-l-J 

F»0.D6996*cXPtRJNR5./7* 1 
RlN-l-KHl 

IFU .cO* KHl) P«0.5*C.lAtXF(RlK*5W7O 
1F<J .EO. KN) P«S 0 P.f( 0 .n 
002 U) •0.2095*1. 013F4C6*P/( KB* TOM 
CA2iJ)«1.7UE*154A.A9fcE*2A* (P-3.t)95E-09> 

2 CONTINUE 
RETURN 
END 


OF 


r K * 

i \r A '• 


o n 


SUBROUTINE HcA.T(TlHE»GBV»PB,r8#OU#KHZ^OEN$»DZ* JH*KN*TZO»BVF>SK*t>l* HEAT 

»SIRJAY*OELUn»i>HO> FEB3 

C0KnuN/BUnLSG/T6Uej#Tfi0P»iei» PRS(7/l9)»OOG HEAT 

»(ia)»oo(iB,i7)»cozsc( id),onz3ab,i7»,cs<iu,iT)*a.aef 17>#ZNU6), heat 

»AZNild>,F0M18),cmi9>i7)^Siril9/17»»C0m9»l7» * HEAT 

COHliON/CUKTIS/PJUlb) ,PHIH(»d)»Cia7)«CZ(lt>I«c<6>i9>lBI HEAT 

«>R (36#17« le > HEAT 

OlHEHStJN CBV(Kn>» PB(JN>KNW TBUHt-KN)# OBIjN^KN)^ 0ENS(KN>, TZO HEAT 

MRN),TPUB>Mi<lU>,CK(2,ie,17) HEAT 

DATA TS/216.7>Z 19. 6# Z2A.Z*229.6> 293.2^2 36.0. 270. 7»2b3. SWAB. Aj> HEAT 

1239. 7 221.0.210.A»201.bW93.2WB6.9.1<l7.lW31.AWSb.2/ HEAT 

JHt - J«-l HEAT 

KNL • KN-l HEAT 

KNP>KN^l 1 KN2-KNA2 HEAT 

J2»J(1L/2 HEAT 

C COMPUTE zonal NEAN TeHPERaTURE (DEVIATION FKCN GLOBAL AVERAGE! HEAT 

DO o5 J-l.JML HEAT 

00 o5 K-l,KNL • HEAT 

65 TB( J»K) •KHZAIPaiJ.K»n*OESS (KfU-PB (J,K )*D£NS (K) >/DZ H^AT 

DO 6b J>1.JHL HbAT 

DO o7 K»l,Ki>L HEAT 

TP(K>-TB(J»K)*TZO(K> HEAT 

67 CONTINUE HEAT 

C CHOOSE APPROPRIATE TROPOSPHERIC TSNPcRATURE HEAT 

TP(KN)-TROP( J) HEAT 

C TP(KN>»TZO(NN) HEAT 

IP(KNP>-TG< J> HEAT 

c cDKPoTt solAR Geometry tactcrs heat 

DAY • Tint/ (29. * 60 . *60. > HEAT 

IDAY-0AYAO.U03 HEAT 

RES - uAY-IOAY heat 

EPS-0. 7*0ELTIM/ (29.*60.*60. > HEAT 

lF(RtS .GT. EPS) GO TO 60 HEAT 

C COMPUfATIOll OF EAHIH-SUH DISTANCE (BERGER, 197»» J. ATMOS. SCI. MEAT 

C 35 PP. 2362-2367) • HEAT 

bC-0. 016722 HEAT 

BETA-SQRT(l.-tC*EC) HEAT 

CM8-(10l.21972+lbO. )*Pl/ieo. HEAT 

kLAHBZ--2.*( (£C/2.*(EC**3 )/8. >A(1.*BETA)ASIN(-0NB) HEAT 

l-£C*tC/9.*(0.5AB£TA)*SIN(-2.A0MB) HEAT 

%* (EC**3)/8.*11. /3.Aa£TA)*SINl-3.*0HB) ) HEAT 

C RLAMBZ IS L0N6IT0DINAL POSITION AT 21 MaRCH (VERNAL ECiOINOX) HEAT 

C ONE YEAR HAS 360 DAYS HEAT 

RLAHBH-RLArtBZ+DAY/360 .*2.*PI HEAT 

RLAf1B-RLA.1BM»(2.PEC-U.25*EC-*3)*SlN(RLAnBrt-OMB) HEAT 

»*(5. /A. )»cC*tCPS IN(2. *RLAMBM-2.*0HB ) HEAT 

!♦ (13. /12. )*EC**3*SiN(3,»RLAMbH-3.A0f1B) HEAT 

V-KLAMB-OrlB HEAT 

RHO • (l.-SC*EC)/(l.*ECACOS(V)) HEAT 

C SEASONAL VARIATION OF SoLAR HEATING (COOLEY AND 80RUCKI (1976) HEAT 

C J. ATMOS. SCI. 33 ?P. 13A7-1356, APPROXIMATION 1) HcAT 

C DELTA-SOLAR DECLINATION HEAT 

DELTA - 0.A093196-SIH(2.*PI*(STR0AY-DAY-U0. )/360. ) HEAT 

DO 30 L-l.JML HEAT 

PHD- (J2*i-L)*10-3 HEAT 

C PHl-IERRESTRIAL LATITUDE HEAT 

PHI - PfiD*PI/l80. HEAT 

TSTAR-IiHE OF jUMitT (CR NEGATIVE TIME OF SUNRISE) HEAT 

ZEK-AVcRAC-u value uF CQS(SZA) BETWEEN -TSTAR AND TSTAR HEAT 

SUN - TAN(OELrA)»IAN(PHI) HEAT 

ISUN - SUN heat 

IF (ISUN) 10,15,20 HEAT 

10 TSTAR - 0 . HEAT 

ZEN - COS (OELTA-PHl ) HEAT 

Co TO 25 ' HEAT 

C HEAT 

15 TSTAR - aZ./PI )*ACOS(-SUN) . HEAT 

ZEN - Si;i(0^rLTA)*SIN(PHI)A(12./(Pl*TSTAR))*CDS(0£LTA)AC0S(PHI)* MEAT 
• SINIPI^TST/VR/IZ.) heat 

GO TO 25 HFAT 

C HEAT 

20 TSTAR - 12. HEAT 

XEN - SIN(OELTA)*SIN(PHI) HEAT 

25 CONTINUE 


he*T 

AZSN • 3>./S0RrU224.*i;eN*ZcN>X.l HEAT 

fOAt • TSTAR/I2. HEAT 

2N(U • ZEN . HEAT 

AZNU) • ATEN HEAT 

FOYIL) • EDAY HEAT 

30 CONTINUE. . HEAT 

C • HEAT 

C COMPUTE TEMEtRATURE-CORRECTEO CURTIS MATRICES FOR EACH ZONE HEAT 

CALI CURHTP/KN,TS^CM> HEAT 

C CONVERT TO JRH VERTICAL COORDINATE LABEL MEAT 

00 6B J^J^lfKNP HEAT 

00 69 lI>i#KNL HEAT 

RU.II» JJ>»CMU»KN2-J J,KN-1 II HEAT 

R (JYJMLfll# JJ)-CN(2»KN2-JJ»KN-1I) HEAT 

69 CONTiNUE HEAT 

r(j»an, jj>-cmi*ks2-jj/Kn> heat 

R ( J^JML>KN>JJ)-CMI2>KN2-JJ»KN) HE AT 

66 CONTINUE HEAT 

66 CONTINUE HEAT 

C COMPUTE GLOBAL RADIATIVE EUUUIBRIUM TEMPERATURE HEAT 

CALL RADEQUITZOfRHO/ JM^KNT HEAT 

IF (TIME.CT.O.) GO TO 7S HEAT 

C COMPUTE STATIC STABILITY FROM TZO HEAT 

N - KNL-1 HEAT 

DO 70 K>2»N HEAT 

70 GBV(K» • aVF*RHZ/<2./7.*TZO(K)/SH*(T2OUFI)-TZ0(K-l n/(2.*DZ) » HEAT 

GBVm • GBVI2) HEAT 

GdVIKNU • CBV(N> HEAT 

WPITE<6»10* HGBV(I) / 1-I.KNL ) HEAT 

75 CONTINUE HEAT 

60 CONTINUE HEAT 

C CCnPuTt HEATING RATE HEAT 

DU 33 L-I/JML HtAT 

00 40 J-1,KNL HEAT 

Tr (J )-I6(L»J)ATZ0(J) HEAT 

40 CONTINUE HEAT 

DO I-1>KNL HEAT 

QBtL,I)-UeLT<TP»L»I»RHO/KN) HEAT 

45 CON riNUE HE AT 

J5 CONTINUE HEAT 

GTIME • I.73E+06 H;AT 

DO 90 K-1#KNL HEAT 

DO 83 fl»JML HEAT 

08 (L/K )-0B( L»K» / (86400. ♦RHZ*DENS(K) )*0Z HEAT 

t*(l.-EXP(-TIME/GTIKE) > HEAT 

35 CONTINUE HEAT 

90 CONTINUE HEAT 

105 FOPHAT (5X/8E15.3) HEAT 

RETURN HFAT 

END heat 



SUBHOUTINE RADCOU(nO>RHO« JH/KN) 

RADCQU 


OlHcNSlON TZ0(17)»A3(Z) 

RAoecy 


DATA lTOT»OT,£PS/2O»0.1#O.l/ 

RAOfcCU 


JrtfJK-l * KNL»RN-X 

R40ECU 


Pl»j.lAli9265A 

R3CECU 


DO 40 HOT 

kAOUU 


> 0. 

RACcCU 


IrtAA IS KNL FOR CLinATOLQOICAL TROPOSPHERIC TEHP£RATURE» 

RAOCCU 


KM FOR RAOIAliVE £OUUlbRIUH TROPOSPHERIC TEMPERATURE 

RAOecU 


IMAK-KML 

RAOECU 


00 I-ItiMAX 

RACtCU 


T20<n • TZ3(n42.*0T 

RAcecu 


DO J*l,2 

raoecu 


T20UJ - T20m-0T 

RACECU 


A3(J) * 0. 

kAOECU 


00 25 L-l.JML 

R AOECU 


THETA • 175-U-i>*10 

kAOECU 


THcTA • THETA»PI/Xeo. 

RACE CU 


W • SIN(TH£TA>*PI/3o. 

RACECU 


A3U)>A3(J)«UELT(TZ0»UURH0>KN)«U 

kAOECU 

25 

CONTINUE 

kACEOU 

30 

CONTINUE 

RAOcCU 


DO • «A3n >-A3(2> )/0T 

RAOECU 


OIFF • A3(2)/DO 

RAOECU 


IN - TZUUT-OXFF 

RACECU 


OTP • OIFF 

RAOEOU 


IF (0TP.GT.2O.) DTP - 20. 

RAOECU 


:F (UTP.LT.-20. ) 0T» • -20. 

RADECU 


TN • T20m-0TP 

PADECU 


UXrr • AaS(OIFF) 

RACECU 


IF (UIFF.GT.SUMJ SUM » OIFF 

PAOECO 


T20(I) - TM 

RACEQU 

35 

CONT INJE 

RAOEOU 


IF (SUK.LT.EPST CO TO 45 

RAOECU 

<r0 

CONTINUE 

RAOECU 


PKINT 65, ITOT 

RAOECU 


STOP X 

RACECU 

4r5 

PKINT 70, IT 

RAOEOU 


WRITE (6, 75) TZO(KN), (TZ0(K),K-1,KNU 

kadeou 


RETURN 

RAOECU 

u 


RAOECU 

65 

FORMAT (5X,*TEMPtRATURE PROFILE FAILED TO CONVERGE AFTER*, 13, 

RACECU 


»* ITERATIONS*) 

RADcCU 

70 

FORMAT (5X,*TEMPERATURE CONVtRCEO AFTER *,I2,* ITERATIONS*) 

RAOECU 

75 

FORMAT (XX,IdF7.2) 

RACECU 


END 

RAOEOU 


FUNCriHH D6LT<TP»L* I>RHO#KfO DEtf 

C0Hrtaii/SUMLtC/TGae»»TPJPIib)>PRS<7»19)#006 0£LT 

l<ie)»00(lB»l7»»0JZSC( lb)»OOZS< ld,i7>»CSae>X7J,CCae»i7)>2N<l8), DELI 

»A2N«18»>FOY.(l«»,C2T(XS» 17)# iZH 19# 17)»C0T( 19, 17) DELT 

COMhDN/CUHTIS/PHI (18) #PHIH( 18 ) #Ca ( 17) # 02 ( 18 ) # £ (6# 19 # 18 ) CPU 

S#Rlio#17# 18) CELT 

C0Mrt0N/02/l)02(17)#CA2 (17) D6LT 

OinENSIOH TP(17), 10 ( 18 ) OcLT 

KMl-KN-l * KNP-KNAl DELT 

8)1L>18 CELT 

C UMb AnENOnENTS# 20 OCTO&En# I960 CELT 

C V-HPAC^WN/KS CELT 

Vli-9b9.96 » V10-J,'<.09.fcA OElT 

Cii-ClO-O, CELT 

C CHOOSt . .•kUPRlAlE TROPOSPHERIC TEMPERATURE DELT 

T0(1)-TG(L) OELT 

T0(2)-TRUP(L) OELT 

C TD(2)-TR(KN) ' OELT 

I>G t, J-3#KNP OELT 

TD(J)-I?( J-2) OfLT 

A CONTINUE OELT 

OU 3 J*1#KNP OELT 

TT*TO(J) pelt 

THIS. (EXP ( V15/300. )-.l. )/(EXP(V15/TT)-l. ) OELT 

THlO- (EXP(V10/300. )-l .)/(EXP(V10/TT)-l. ) 05 LT 

C15*C15aK)L^JHL#1#J )ATH15 CELT 

C10-C10+R(L# I#J)*TH10 OELT 

3 COhTlNUc DELT 

CC02»C13fC10 OELT 

COT (L# I ) »CC02 * OELT 

CTS*-=CS(L#1 >*cXP(-15oO./TP( I) ) OELT 

CTG-tG(L, n*(EXP(^1500,/TG(L))-EXP(-1500. /TP(1) ) ) DELT 

COZ0NE»CTS»CTC DELT 

CZT (L,I ) • COZONE OELT 

OELT-CCU2ACOZONE OELT 

S0»0. CELT 

SZT(L>I)»SO DELT 

ZEN-ZN(v) DELT 

A2£N«AZN(L) OELT 

POAY-FDr(L) OELT 

1F(P0AY .LT. 0.0001) 60 TO 25 DELT 

DZ • 5000. OELT 

SH • 7000. DELT 

Z • (1-1)*DZ DELT 

AlKDEN»1.0136»0b*0.0o996AEXP(-Z/Srt)/(2.87E*06*TP( I) ) OELT 

1F(I .lO. KN) AlROEN-1. 103EY06A0. 3162/(2. 87E«06AT0(2) ) DELT 

TAU • 10. ■ CELT 

CLOUD • 0.996 OEtT 

RG • 0.3 DELT 

ICLOUO-KN * DELT 

C HEATING BY THE A8SOKPIION OF SOLAR RAOIATiON BY OZONE OELT 

C FKON LACIS A^0 HANSEN# J ATNOS SCI 31 118-133 OELT 

C CLEAR SKY DELT 

RAB > 0. 219/(1. *0.U16*ZEN) OELT 

hOB • W.199 OELT 

Rul - RAB*(1.-RAB)A|1.-R0B)ARC/(1.-R0BARG) OELT 

KB - RBI OELT 

C CLOUDY SKY OELT 

HAB - 0.13*TAU/ (1.+0.13ATAU) OELT 

ROB ■ RAQ ‘ OELT 

RBI - RAQf (l.-RAd)*(l.-ROfl)ARG/(l.-ROB*RC) CELT 

UT - AZEN*OUZS(L#ICLOUO)^1.9*(OOZSIL#ICLOUO)-QOZS(L#X)) DELT 

A1 • OZUV(UT) OELT 

UT ■ AZEn»OOZS6(L)«1.9*(OCZS/«(L)-OOZS(L#I) > DELT 

A2 • OZUV(UT) CELT 

U - AZEN*OOZSIL»l) DELT 

A3 • OZUV(U) OELT 

CP-l.tA07 A SOLFLX-0. 1365E*07 DbLT 

SC3-FDAYA86900.AS0LFLX*ZEN*Q0)L#i) /(CP*AIR0EN)*(AZENAA3* CELT 

SCL0U0*1.9*RB1AA1A(1.-CL0U0)A1.9ARB*A2) DELT 


(>K N AJ. 1 A*'-? 

OF P(K)K 


U-A..eN*CA2U» DtLT 

^.02UV«U> DEI.T 

Swi2-A»DO2(i)*8!>«OO.»FDA¥/<AI«0EN*CP) OELT 

i0«tS03»S02)/(RH0*RH0) OEtT 

Smt,n»SO Ot-LT 

25 DtLT»OtLT*SO OELT 

RETURN . OELT 

ENO OELT 


(Changes in main program) 

RCA0(2*19)M (( n* 

143 F0gHAT(3(SE20. ;3/)»AE20.13) 

ReA0(2*i.<iMKAPn>» ir-l»KN> 

READI2»2U0 I T2 0(KN»« (TZO(li« !>l*KNO 
191 FORHATtFT.A) 

READ (2>ia0» ( (PRS<UJ)> 

REA0(2«19t>) QOC»<-iO(L»RN»»L»l. JHL». H30a»K)»L-l,JKLr»K-l*KNL) 
REA0<2>199) aozsc* «»J02 S KN) » L»l» JHU » ( ( OOZ S ( L/K ) / L -I , JHL )» 
*K-1*KNL> 

RLA0<2»195MCISU#KS>.l*l» jnL>»«(CTSa»K>»L-l»JM(,»,K-l»KNl) 
READ(2/19»MCG (L, KMI > L - 1* JHU » (ICC ( If K t Ir JHL ) / K- 1,)^NL I , 

ISO FOkRAT (AF13.A) 

195 FORRAT (9cl0.3> 

200 FORRAT (9F),0.3> 


DO 20 K>1>KN 


if (KoGTsKNt) <50 TO 20 

C NEWTOMUH COOLING CO£FfICIEMTS fOR EDDIES 
KAP<K)-KAP<K>* 0 T/ 86 ^UJ. 

20 CONTINUE 


IRAQ • 0 
1RCT»12 


C COMPUTE ZONAL MtAH TENPERATURE (OEVIATION fROH GLOdAL AVERAGE) 
00 6^ J-1#JHL 
00 6A K-lfKNL 

6^ TB( JiK) -RHZMPdi J#K^l)^OENS(Kfl)-PBU#K)AOENS<K) ) / OZ 
OU 66 J-1#JHL 
00 67 K-X/KNL 
TP(K)*TB(J>K)f TZO(K> 

67 CONTiriUE 

C CHOOSE APPROPRIATE TRQPOSPHcRlC TEMPERATURE 
TP(KN)-TROP< J) 

C TP(KN)«TZU(KN) 

TP(KNP)-TG(J) 

C ClJMPUTc TErtPcRATURE-CORRECTEO CURTIS MATRICES FOR EACH ZONE 
CALL CURT! IP^KN, TS#CMAT) 

C CONVERT TO JRH VcRTICaL COOKOINATE LAPEL 
00 68 JJ-l/KNP 
DO 69 1I-I#KML 

CURMAT (J, X r# JJ )-CMAT( 1 # KK‘2- J J# KN-I I ) 

CURMATI J^JML^i I, J J ) -CMAT (2 / KN2-J J # KN- II I 
69 CONTINUE 

CURMAT IJ/KN# JJ )-CMAT( l#KN2-JJ^Kin 
i.uRhAT( JiJML#KN/ J J ) -CMAT(2#KK2--J J^KNI 

68 CONTINUE 
66 CONTINUE 

lao format (6F15.AI 
195 FORMAT ^9E10.3) 

200 FORMAT (9F10.3) 


WAVE2 

WAVE! 

WAVE2 

WAVF2 

NAVF2 

WAVE2 

UAVE2 

WAVE2 

WAVI2 

WAVE2 

WAVEI 

WAVF2 

WAVE2 

WAVE2 


WAVE2 


yAVE 2 

VAVE2 

WAVE2 

NAVE2 


WAVE2 

FCB3 


WAVE2 

WAVE2 

WAVE2 

WAVF2 

WAVE2 

WAVE2 

WAVE2 

WAVE2 

WAVE2 

WAVE2 

VIVE 2 

UAVE2 

WAVE2 

WAVE2 

WAVE2 

WAVE2 

WAVE2 

WAVE2 

VAVE2 

WAVE! 

WAVE2 

WAVF2 

VAVE2 

WAVE2 

WAVE2 

WAVE2 

W4VE2 


C COHPUTt INITIAL TROPOSPHERIC AMO GROUND TCNPCRATURES WAVFT 

J2•JNl/^ WAVC2 

DO 96 L-1»J2 NAVE2 

PMO-L*XU-S UAVEZ 

TANP«(0. A7^*PHO-A. M2. ¥*VC2 

ARC-2.*P1/36U.* (STROAY-leb.l NAVP2 

IC0X»54.O5*(CO5(PHD*PI/iBO. )-0.22I-lS. WAVE! 

TGIJ2*X-U-UO)lATA(1P»C0S (ARG> WAVE 2 

TC< J2*X-l ) "TCI J2A1-L M27 J. 16 WAVE2 

TC<L4J2( -TEQJ<-TANP*C05 ( ARC> WAVE2 

TGUfJEl-TCa* J2>*273.16 WAVI2 

TaMP-^*.3«»XA1PH0-o.M2. WAVI2 

IE<PHO .CT. AO. I TA.IP.6. WAVE2 

ARG*2.*?l/360.A(STROAr-lOO.I WAVE! 

TE0X«2b.X3*(C0S<PHn»PI /l6u.)-0.26)-5A. WAVE! 

l,E(PHO .LT. 30. » TEQX*-32. WAVE2 

rPOP( J2AI-U -rEOXATAHP*COS( ARC» HAVE2 

TPUP( J2*i-L» -TROPl J2M-LM273. 16 WAVE2 

TROP< J24L) •TtOX-TAHP*:0S <ARC» WAVE2 

lRaP<J2*U«TR0PU2AL)*27 3.16 VAVC2 

96 COrtriHUE WAVE2 


